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About  Tms  Report 

This  report  is  one  of  ten  individual  studies  conducted  for  the  Flathead  Basin  Forest  Practices/ 
Water  Quality  and  Fisheries  Cooperative  Program.  The  Cooperative  Program  was  administered  by 
a  Coordinating  Team  representing  the  Montana  Department  of  State  Lands  Forestry  Division,  the 
Flathead  National  Forest,  Plum  Creek  Timber  Company,  L.P.,  the  Montana  Department  of  Fish, 
Wildlife  and  Parks,  the  Montana  Department  of  Health  and  Environmental  Sciences'  Water  Quality 
Bureau,  the  University  of  Montana,  and  the  Flathead  Basin  Commission. 

The  Cooperative  Program's  specific  objectives  were  (1)  to  document,  evaluate,  and  monitor 
whether  forest  practices  affect  water  quality  and  fisheries  within  the  Flathead  Basin,  and  (2)  if 
detrimental  impacts  exist,  to  establish  a  process  to  utilize  this  information  to  develop  criteria  and 
administrative  procedures  for  protecting  water  quality  and  fisheries. 

The  ten  individual  studies  included  the  evaluation  of:  (1)  specific  practices  at  the  site  level,  (2) 
accumulation  of  practices  at  the  watershed  level,  (3)  general  stream  conditions,  (4)  water  quality 
variables  relative  to  levels  of  management  activity  in  small  watersheds,  (5)  fish  habitat  and 
abundance  relative  to  stream  variables  influenced  by  forest  practices  at  the  watershed  level,  (6)  long- 
term  changes  in  large- stream  dynamics  related  to  historical  records  of  natural  and  man-related 
disturbances,  and  (7)  changes  in  lake  sediments  relative  to  historical  records  of  natural  and  man- 
related  disturbances.  A  Final  Report  was  developed  which  contains  summaries  of  each  of  the  studies, 
a  set  of  summary  conclusions  and  recommendations,  and  a  formal  response  to  the  recommendations 
by  the  land  management  organizations  which  administered  the  Cooperative  Program. 
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EXECUTIVE  SUMMARY 

Geographic  Information  Systems  (GIS)  are  rapidly  gaining 
recognition  as  powerful  tools  for  land  use  planning  and 
management.   This  has  been  accompanied  by  an  equally  rapid 
development  of  GIS  Applications  -  methodologies  developed  to  make 
use  of  the  GIS  capability  to  link  tremendous  quantities  of 
spatial  and  nonspatial  data. 

This  report  documents  the  development  of  one  of  those 
applications.   Funding  for  this  study  was  not  obtained  from  the 
Flathead  Cooperative,  but  rather  from  a  U.S.  Department  of 
Agriculture,  Mclntire-Stennis  grant.   The  majority  of  development 
work  was  done  on  the  University  of  Montana,  School  of  Forestry's 
PAMAP  GIS.   However,  the  application  was  developed  to  be 
independent  of  any  particular  software  package.   We  worked  in 
close  communication  with  the  Flathead  National  Forest  MURIS  staff 
to  make  sure  that  the  "layers"  of  information  used  in  the 
application  were  either  available  or  readily  synthesized  from 
existing  information  on  their  system. 

The  body  of  this  report  is  the  Master  of  Science  thesis 
prepared  and  successfully  defended  by  Kenneth  J.  Lull. 
Refinements  on  the  application  continued  after  the  defense  of  his 
thesis,  however.   Appended  to  the  report  is  the  completed  Risk 
Matrix  which  is  central  to  the  application,  and  a  paper  that  was 
presented  at  a  national  symposium  on  computer  models  and 
applications.   The  paper  outlines  the  methods  for  assessing 
potential  impacts  and  developing  thresholds.   Similar 
applications  have  been  developed  for  use  on  agricultural  lands. 


but  to  our  knowledge,  this  is  the  first  to  be  developed  for  use 
on  forested  lands. 

The  application  has  yet  to  be  tested  on  the  Flathead 
National  Forest  MURIS,  but  that  will  be  accomplished  soon. 
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Lull,  Kenneth  J.,  M.S.  June.  1990 

A  GIS  application  for  assessment  of  nonpoint  source 
pollution  risk  on  managed  forest  lands.  92  pp. 

Director:  Dr.  Donald  F.  Potts 

Legislation  In  recent  years  has  required  states  and 
federal  agencies  to  control,  assess,  and  monitor  nonpoint 
source  pollution  resulting  from  land  management.   A 
methodology  to  assess  the  risk  of  water  quality  degradation 
resulting  from  timber  harvest  activities  using  geographic 
Information  systems  Is  presented. 

An  erosion-impact  matrix  Is  developed  assigning  risk 
values  to  differing  s I  I v I cul tura I  practices  based  on 
land-type  Information.   Data  assimilation  and  spatial 
analysis  for  the  study  area  is  conducted  using  GIS.   The 
risk  matrix  Is  developed  using  a  modified  Delphi  technique. 
Results  from  two  pilot  studies  in  western  Montana  suggest 
the  utility  of  GIS  In  watershed  management.   The  capability 
of  advanced  GIS  technology  allows  editing,  storing,  and 
spatial  analysis  to  be  conducted  In  a  timely  and  cost 
efficient  manner. 
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ORGANIZATION 

There  are  two  major  portions  of  this  study,  they  will 
be  discussed  Independently. 

Chapter  I  -  Erosion  Hazard  Matrix  Development  - 
includes  detailed  explanation  of  erosion  hazard  matrix 
construction,  definition  of  matrix  components,  and  product 
risk  values  for  forested  watersheds  In  western  Montana. 

Chapter  II  -  Geographic  information  System  (GIS) 
risic  assessment  application  -  describes  the  process  of 
combining  spatial  and  nonspatlal  data  simultaneously  in 
conjunction  with  erosion  risic  values.   Includes  description 
and  results  for  two  pilot  watersheds  in  western  Montana. 


CHAPTER  1 
EROS  I  ON- IMPACT  MATRIX  DEVELOPMENT 


INTRODUCTION: 

Cumulative  watershed  effects  can  be  thought  of  as  the 
"total"  impact  on  aquatic  resources  resulting  from  a 
combination  of  varying  land-use  activities  upstream.  A 
primary  pollutant  from  these  nonpoint  source  activities  is 
sediment.  Nonpoint  source  (NPS)  pollution  generally  results 
from  diffuse,  uncontrolled  surface  runoff  caused  by  various 
activities.  Including  agriculture  and  silviculture 
operations  (Myers  &  Wise,  1989).   in  a  1986  Environmental 
Protection  Agency  report  to  congress,  NPS  pollution  was 
responsible  for  failure  to  meet  water  quality  standards  in 
76%  of  the  iaice  acres  reported  by  states,  65%  of  the  stream 
miles,  and  45%  of  the  estuarine  waters. 


Page  2 

The  subtle  nature  of  NPS  pollution  presents  significant 
difficulties  to  agencies  charged  with  monitoring  and 
management  (Myer  &Wise,  1989).   impacts  from  Individual 
disturbance  activities  may  be  minimal,  but  collectively  and 
cumulatively  significant.   Hence,  NPS  pollution  is  among  the 
nation's  most  serious  natural  resource  problems  (Myers, 
1986). 

Forest  management  and  multiple-use  encroachment  Into 
forested  watersheds  may  result  in  negative  Impacts  to  both 
soil  and  water  resources.   For  example,  construction  of  a 
logging  road  compacts  bare  soil,  which  reduces  infiltration 
and  percolation,  and  can  change  runoff  processes. 
Harvesting  and  site  preparation  activities  also  expose  bare 
mineral  soli  which  In  turn  can  lead  to  soil  erosion  and 
stream  sedimentation  during  heavy  rainstorms.   If  these 
events  combine  with  increased  peak  flows  because  of  timber 
harvest  (which  reduces  evapotransplrat Ion)  during  a  heavy 
rain  or  snowmelt  event,  mass  wasting,  channel  erosion  and 
aggradation  can  accelerate  beyond  natural  levels  (Coburn, 
1989).   Potential  Impacts  Include  a  decline  In  fish  spawning 
habitat.  Increased  costs  associated  with  watershed 
rehabilitation,  and  financial  burdens  due  to  flood  damage  on 
public  and  private  property. 
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As  a  result  of  public  and  scientific  concern  over  water 
quality  protection,  various  federal  laws  have  been  enacted. 
An  Important  precursor  of  the  major  legislation  of  the 
1970's  was  the  Multiple  Use  and  Sustained  Yield  Act  of  1960 
(Wilkinson  &  Anderson,  1985).  The  guiding  principle  of 
sustained  yield  was  defined  as  "the  achievement  and 
maintenance  in  perpetuity  of  high-level  annual  yield  or 
regular  periodic  output  of  various  renewable  resources  of 
the  national  forests  without  Impairment  of  the  productivity 
of  the  land"  (Coburn,  1989). 

The  National  Environmental  Policy  Act  of  1969  (NEPA- 
P.L.  91-190),  requires  federal  agencies  to  assess  the 
cumulative  and  long  term  effects  of  proposed  land-use 
actions.  Additionally,  Section  102  (2c)  calls  for  the 
preparation  of  an  environmental  Impact  statement  on  projects 
that  may  significantly  affect  the  quality  of  the  environment 
as  part  of  the  planning  process. 

The  next  significant  piece  of  federal  legislation  to 
set  requirements  for  consideration  of  cumulative  effects  was 
the  1977  Clean  Water  Act  Amendments  (P.L.  95-217),  This 
act  established.  In  Section  208,  a  process  for  controlling 
nonpoint  sources  of  pollution  and  required  states  to  prepare 
water  quality  standards  that  are  determined  by  the  highest 
"beneficial  uses"  of  the  water  In  question  (Coburn,  1989). 
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Section  319  of  the  Clean  Water  Act  of  1987  (P.L.  100-4) 
additionally  requires  Individual  states  to,  as  a  minimum, 
initiate  plans  that  encompass  three  primary  objectives:  1) 
assessment.  Identification  ,  and  definition  of  water  quality 
problem  areas  and  their  causes  -  point  and  nonpoint;  2) 
targeting  action  priorities  applied  to  assessed  waters  based 
on  indicators  such  as  public  health,  environmental  risks, 
and  value  of  aquatic  habitat;  and  3)  Identification  and 
Implementation  of  Best  Management  Practices  (BMPs)  for 
rehabilitation  and  prevention  of  further  degradation  of 
targeted  waters  (Meyers  &  Wise,  1989). 

In  an  effort  to  curb  the  detrimental  Impacts  of  NPS 
pcilutlon  and  comply  with  federal  legislation,  cumulative 
watershed  effects  analysis  Is  necessary.  Cumulative 
watershed  effects  (CWE)  analysis  is  basically  an  advanced 
means  of  controlling  nonpoint  source  pollution 
(Coburn,1989) .  The  process  Includes  predicting  Impacts  that 
may  deteriorate  water  quality  which  might  be  missed  if 
planning  were  carried  out  only  at  the  project  proposal 
level.  By  recognizing  that  a  watershed  Is  a  fluvial  system 
and  an  ecological  entity,  cumulative  watershed  effects 
analysis  assumes  that  hydrologic  effects,  erosion  processes, 
and  biological  responses  are  considered  from  the  outset 
(Coburn,  1989).  Long  term  benefits  Include  maintenance  of 
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water  quality  and  cold  water  fisheries. 

Definition  of  tine  term  CWE  has  been  slow  and  Is  still 
not  universally  accepted.  Cumulative  effects  have  been 
called  the  UFOs  of  hydrology  (Grant,  1985).   Hasklns  (1986) 
described  CWE  as  the  additive  and/or  synergistic  effects  of 
land  management  activities  on  water  quality  and  beneficial 
uses,  which  occur  away  from  the  site  of  primary  development, 
and  are  transmitted  to  the  fluvial  system.  The  current 
confusion  about  cumulative  effects  stems  In  a  large  part 
from  a  lack  of  conceptual  models  and  research  tools  for 
analyzing  the  complex  nature  of  drainage  basin  response  to 
disturbance.  Without  such  models  or  techniques,  the  term 
cumulative  effects  Is  sometimes  Invoked  to  explain 
phenomena  we  observe  but  are  unable  to  attribute  to  a 
specific  cause  (Grant,  1985). 

Forest  land  managers  are  concerned  about  the  potential 
cumulative  effects  of  multiple  management  activities  over 
time  and  space  within  a  watershed  on  the  downstream  aquatic 
ecosystem  (Klock,  1985).  The  most  appropriate  approach  to 
determine  the  cumulative  effects  of  forest  practices  on 
water  quality  would  be  a  large  watershed  study.  However,  In 
light  of  such  vast  financial  and  resource  commitments, 
watershed  analysis  at  this  scale  Is  often  not  very 
practical.  Thus,  an  alternative  approach  is  to  develop  a 
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"watershed  cumulative  effects  analysis  model  which  best 
reflects  the  multitude  of  potential  downstream  Impacts 
forest  practices  may  generate"   (Klock,  1985).   Rlckert  et 
a!  (1978)  suggests  the  use  of  an  erosion-hazard,  or  risk 
analysis  model.  In  order  to  predict  potential  NPS  pollution 
problem  areas. 

Coburn  (1989),  recommends  a  computerized  database, 
using  state-of-the-art  geographic  information  systems  to 
enable  long-term  tracking  and  maintenance  of  sensitivity  and 
land  disturbance  Inventories. 

A  geographic  Information  system  is  designed  to  accept 
large  quantities  of  spatial  data  derived  from  a  variety  of 
sources.   Such  a  system  Is  also  designed  to  store,  retrieve, 
manipulate,  analyze,  and  display  those  data  according  to 
user-defined  specifications  (Marble  and  Penquet,  1983). 
Geographic  Information  systems  (GIS)  are  becoming 
Increasingly  important  in  natural  resource  management 
(Morgan  and  Neiepa,  1982).  Using  a  geographically 
referenced  system  with  soil  erosion  and  non-point  source 
pollution  models  can  account  for  every  parcel  of  land  and, 
therefore,  allow  researchers  to  evaluate  large  geographic 
areas  in  much  less  time  (Pel  I  1 1 ier , 1985) . 
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OBJECTIVES 

In  keeping  with  the  ever-expanding  needs  of  natural 
resource  management  and  the  advanced  capabilities  of 
computer  graphic  and  database  management  systems,  the  goal 
of  this  project  was  to  develop  a  nonpoint  source  risk 
analysis  procedure  for  use  in  geographic  Information 
systems.   The  following  methodology  Is  based  on  a  nonpoint 
source  risk  analysis  procedure  developed  In  Oregon  (Rlckert 
et  al.  1978:  Brown  ill  et  al.  1979),  but  differs  in  that  all 
cartography  and  spatial  analysis  will  be  done  through  the 
use  of  GIS.  There  are  two  major  components  to  this  study: 

1)  the  development  of  an  erosional-lmpact  risk  matrix;  and 

2)  the  Implementation  of  the  risk  procedure  using  GIS  on  two 
pilot  study  areas  In  western  Montana. 

MATRIX  CONSTRUCTION 

Maps  and  associated  topographic  attributes  have  long 
been  used  by  resource  managers  in  decision  making  and 
land-use  planning.  The  utility  of  the  land-type  map  can  be 
greatly  enhanced,  however,  with  the  use  of  a  scheme  for 
ranking  the  erodlblllty  of  different  land  surfaces  under 
different  management  scenarios  for  a  given  watershed.  The 
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scheme.  In  this  case,  Is  an  erosion- Impact  matrix  which 
provides  a  systematic  basis  for  malcing  estimates  of  the 
relative  Impact  or  risk  of  NPS  water  quality  degradation  of 
human  activity  on  different  types  of  terrain.   In  the  matrix 
developed  by  Brown  ill  et  al.  (1979),  the  horizontal  axis  is 
composed  of  order-of-magn I  tude  NPS  rlsi<  factors  for  geology 
and  slope  combinations.   The  vertical  axis  Is  composed  of 
order-of-magn I tude  erosionai-impact  factors  for  selected 
land-use  activities. 

The  erosion-Impact  matrix  developed  for  this  study  is 
similar  in  that  the  land-use  activities  are  on  the 
horizontal  axis  and  slope  and  geology  are  on  the  vertical 
axis  (See  Figure  1).  The  matrix  was  modified  to  reflect 
user  needs  based  upon  geography  and  resource  concerns  more 
specific  to  western  Montana,   it  Is  suggested  that  matrix 
construction  be  based  on  the  socio-economic  importance  of 
the  basin  for  a  given  region  (Brown  III  et  al.  1979). 
Therefore,  It  was  decided  that  the  emphasis  of  the  erosion 
-Impact  matrix  for  this  study  focus  on  the  Implications  of 
forest  harvest  activities.   The  same  matrix  is  used  for  both 
pi  lot  watersheds. 

The  use  of  the  eros lona l-lmpact  matrix  allows  the 
generation  of  a  new  family  of  map  products  by  the  GiS, 
derived  from  the  original  slope/substrate  maps  and  based 
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FIGURE  1:  MONTANA  EROSION-IMPACT  MATRIX 
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upon  modification  of  land-use  activities.  A  nonpoint  source 
risk  map  founded  on  existing  watershed  conditions  may  also 
be  developed,  given  some  basic  assumptions  about  recovery 
periods  following  disturbance.  Both  of  these  applications 
will  be  discussed  in  the  following  sections  and  in  Chapter 
2. 

A  major  problem  with  the  use  of  an  erosion  matrix 
centers  on  the  order-of-magn I tude  values  within  the  body  of 
the  matrix  Itself.   It  Is  relatively  easy  to  obtain 
agreement  on  the  overall  parameters  which  comprise  the 
horizontal  and  vertical  administrative  data  (I.e. 
slope/geology  and  land-use  activities),  however,  placing 
values  Internally  within  the  matrix  Is  extremely  difficult. 
In  an  attempt  to  solve  this  problem.  It  was  decided  that 
"best  professional  Judgement"  would  be  used  and  that  an 
approach  similar  to  the  Delphi  Technique  be  applied. 

THE  DELPHI  TECHNIQUE 

As  previously  described,  the  Oregon  matrix  for 
management  practices  and  land-type  erosion  factors  was 
modified  for  use  in  western  Montana.   Risk  potential  ratings 
for  the  Montana  eros iona i-lmpact  matrix  were  obtained 
through  the  Delphi  method.  This  was  done  in  coordination 
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with  the  Montana  Riparian  Association. 

A  primary  objective  of  a  Delphi  study  Is  to  obtain 
consensus  about  Judgmental  Information  from  a  group  of 
expert  respondents  (DaJanI  et  al.,  1979).  Typically, 
members  of  a  Delphi  panel  are  knowledgeable  In  the  special 
field  of  application  and  they  try  In  a  discuss  Ion-type 
procedure  to  contribute  their  experience  In  order  to  solve 
problems  or  answer  certain  questions  (Rauch  1979,  Khan 
1989). 

A  similar  procedure  was  used  In  the  Fall  of  1988  to 
obtain  eroslonal-lmpact  matrix  values  for  the  study  area. 
Experts  In  soil  science,  hydrology,  silviculture,  ecology, 
and  fisheries  were  Included  In  the  Delphi  panel.  Each  panel 
member  was  sent  an  Initial  survey  and  Information  paclcet 
that  explained  Delphi  procedures,  matrix  construction,  and 
the  objectives  of  the  study.  A  blanic  eroslonal-impact 
matrix  was  included  In  each  mailing  and  the  panel  members 
were  asiced  to  evaluate  variable  land-use  activities  with 
respect  to  given  slope  and  parent  material  soil  erodlblllty 
categories  In  terms  of  the  potential  rlsl<  of  NPS  pollution 
production  and  cumulative  watershed  effects.  Each  activity 
was  given  a  rating  from  1  to  5  with  the  following  risk 
values  assigned:  a)  5-  very  high  risk;  b)  4-  high  risk;  c) 
3-  moderate  risk;  d)  2-  low  risk;  and  e)  1-  very  low  risk. 
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For  example,  a  cable-logging  seed  tree  harvest  with 
full  suspension  on  a  soils  derived  from  alluvial  parent 
materials  might  receive  a  risk  erosion-Impact  rating  of  3.0, 
moderate.   The  objective  In  using  the  Delphi  format  Is  to 
try  to  obtain  consensus  from  all  panel  members  that  this 
activity  does  indeed  pose  a  risk  value  of  3  In  terms  of 
potential  NPS  pollution  production.  Once  all  cells  of  the 
matrix  have  been  filled  In,  they  can  be  used  to  evaluate  the 
potential  risk  of  future  activities,  as  well  as  those 
already  existing,  on  an  areal  basis. 

There  are  several  very  Important  factors  which  must  be 
evaluated  before  a  Delphi-like  approach  can  be  effectively 
utilized.   First,  and  probably  the  most  important,  is 
anonymity.  Anonymity  Is  essential  to  guarantee  that  Ideas, 
concepts  and/or  arguments  are  not  influenced  by  other  panel 
members  (Rauch,  1979).  The  director  of  the  Delphi  study 
should  Insure,  at  least  in  theory,  that  the  results  lead  to 
an  unbiased  and  complete  picture  and  solution  to  the  problem 
at  hand. 

In  general  terms,  the  Delphi  procedure  exists  In  2 
distinct  forms.   The  most  common  version  Is  commonly 
referred  to  as  a  "Delphi  Exercise".   This  typically  requires 
a  "paper-and-penci I"  drill  on  behalf  of  the  moderator.   In  a 
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Delphi  exercise,  a  small  monitor  team  designs  a 
questionnaire  which  is  sent  to  a  larger  respondent  group. 
The  respondent  group  Is  usually  given  at  least  one 
opportunity  to  reevaluate  Its  original  answers  based  upon 
examination  of  the  group  response.  This  form  of  Delphi  is  a 
combination  of  a  polling  procedure  and  a  conference 
procedure  which  attempts  to  shift  a  significant  portion  of 
the  effort  needed  for  individuals  to  communicate  from  a 
larger  respondent  group  to  a  smaller  monitor  team.   This  Is 
essentially  the  procedure  that  was  used  to  obtain  erosion- 
impact  risk  values  for  our  matrix.   However,  interspersed 
between  mailings  was  a  meeting  of  all  panel  members  to 
discuss  the  previous  rounds  results. 

The  second  and  newer  form  Is  often  called  "Delphi 
Conference".   It  replaces  the  monitor  team  with  a  computer 
which  has  been  programmed  to  carry  out  the  compilation  of 
group  results.  This  latter  approach  has  the  advantage  of 
eliminating  the  delay  caused  In  summarizing  each  round  of 
Delphi,  thereby  turning  the  process  Into  a  real  time 
communications  system  (Linstone  and  Turoff,  1975). 
Discussions  on  computer  applications  of  Delphi  can  be  found 
within  articles  published  elsewhere,  namely  Johansen  et  al. 
(1974),  Johansen  and  Shuyier  (1975),  and  Turoff  (1972).   The 
University  of  Montana  does  not  have  the  automated  capability 
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or  the  software  to  conduct  computer  Delphi  procedures  and 
therefore  the  procedure  used  In  this  study  Is  classified  as 
a  "Delphi  Exercise".   Further  elaboration  on  computer  Delphi 
methods  Is  beyond  the  scope  of  this  text. 

Delphi  procedures  usually  consist  of  four  major  phases. 
The  first  phase  Is  characterized  by  exploration  of  the 
subject  under  discussion,  wherein  individual  members  of  the 
panel  and  monitor  team  contributes  additional  Information 
that  he  of  she  feels  pertinent  to  the  Issue.   The  second 
phase  involves  the  process  of  reaching  an  understanding  of 
how  group  members  view  the  issue  (i.e.  were  the  members  in 
agreement  or  disagreement  on  what  they  mean  by  relative 
terms  such  as  Importance,  desirability,  or  feasibility),   if 
there  is  significant  disagreement,  than  It  must  be  explored 
In  the  third  phase  to  bring  out  the  underlying  reasons  for 
the  differences  and  possibly  to  evaluate  them.   The  last 
phase,  a  final  evaluation,  occurs  when  all  previously 
gathered  information  has  been  Initially  analyzed  and  the 
evaluations  have  been  fed  bacic  for  consideration  (tinstone 
and  Turoff ,  1975). 

Creating  the  panel  of  experts  is  usually  the  first 
task.   It  Is  prudent  to  choose  qualified  but  diverse  panel 
members  In  order  to  cover  all  aspects  of  the  problem  at 
hand.   For  creating  a  successful  mix  of  panelists,  Scheeie 
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(1975)  suggests  using  three  types  of  personalities.  These 
include:  1)  stake-holder,  those  who  are  or  will  be  directly 
effected;  2)  experts,  those  who  have  an  applicable 
speciality  or  relevant  experiences;  and  3)  facilitators, 
those  who  have  skills  in  clarifying,  organizing, 
synthesizing,  st imulat ing. . . .p lus.  Individuals  who  can 
supply  alternative  global  views  of  the  culture  and  society. 

Schiebe  et  al.  (1975)  describes  a  possible  step  by  step 
procedure  which  In  theory  can  be  applied  to  most  Delphi  case 
studies.  The  process  follows: 


Step  1:  The  process  administrator  (PA)  explains  the  Delphi 

method  to  the  panel. 
Step  2:  PA  administers  a  Pre-Delphi  questionnaire.   This 

helps  to  draw  a  personnel  profile  of  panel  members 

and  may  be  used  to  facilitate  the  removal  or 

addition  of  panel  members. 
Step  3:  Panel  fills  out  the  questionnaire. 
Step  4:  The  PA  describes  a  hypothetical  situation  to  panel. 
Step  5:  Panel  defines  objectives. 
Step  6:  PA  categorizes  objectives  and  panel  assigns  ratings 

or  ranks  the  goals. 
Step  7:  The  PA  compiles  the  1st  round  rankings  and  returns 

feedback  to  the  panel. 
Step  8:  The  panel  rates  goals  for  round  two  with  written 

comments. 
Step  9:  The  PA  compiles  ratings  and  written  comments  and 
returns  all  information  to  the  panel. 


This  procedure  repeats  Itself  up  to  four  consecutive  but 
distinct  rounds.   Between  each  round  it  is  suggested  that 
the  panel  have  a  formal  meeting  under  the  direction  of  the 
PA  to  discuss  the  previous  rounds  results.  This  leads  to 
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group  Interaction  and  allows  panel  members  to  share  personal 
experience.  At  the  end  of  the  fourth  round,  the  PA  complies 
and  analyzes  all  ratings  for  Individual  rounds.  Correlation 
of  the  questionnaires  and  publication  of  results  Is  now  In 
order.  A  graphical  representation  of  this  process  as 
published  by  Schlebe  et  al.  (1975)  can  be  seen  in  Figure  2. 

Termination  of  Delphi  processes  is  often  not  clearly 
defined  In  the  literature.  Dejani  et  al.  (1979)  presents  a 
hierarchical  stopping  criterion  for  Delphi  studies. 
Assuming  that  respondents  have  exhibited  stability  in  their 
responses  In  two  consecutive  rounds,  termination  of  the 
Delphi  procedure  may  take  place,  but  only  If  stability 
manifests  In  any  of  the  following  ways: 

1.  Consensus:  occurs  when  unanimity  is  achieved 
concerning  any  given  issue.  When  consensus  results  the 
study  may  be  terminated. 

2.  Majority:  occurs  when  more  than  50%  of  the 
respondents  exhibit  consistency.   When  a  majority  occurs  and 
is  coupled  with  an  apparent  agreement  among  the  minority 
respondents,  that  the  study  may  be  terminated. 

3.  Bipolarity:  occurs  when  respondents  are  equally 
divided  over  an  Issue.  When  bipolarity  occurs  one  should 
determine  the  nature  of  the  stability  among  the  two  bipolar 
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groups.  A  decision  must  then  be  made  as  to  whether  to 
terminate  or  rewrite  the  particular  question.   If  the  latter 
choice  Is  made,  a  new  round  Is  administered  using  the 
rephrased  question.  This  question.  In  turn,  must  pass  the 
stability  test  before  it  is  dropped  from  the  study. 

4.  Plurality:  occurs  when  the  larger  portion  of  the 
respondents  (but  less  than  50%)  reach  agreement.  When 
plurality  occurs  one  should  checic  for  consensus  within  the 
pluralistic  group  and  for  the  nature  of  stability.  If  any, 
among  other  Individual  respondents  or  clusters  of 
respondents.   If  stability  Is  not  established,  a  new  round 
of  questions  is  administered.   If,  on  the  other  hand, 
stability  is  established,  there  Is  a  choice  of  terminating 
the  particular  question  or  of  rephrasing  It  and  including  it 
In  the  following  rounds. 

5.  Disagreement:  occurs  when  each  respondent  maintains 
views  independent  of  each  other  respondent,  such  that  the 
responses  cannot  be  brought  Into  consonance.  Whenever 
stable  disagreement  is  achieved  for  a  given  question,  the 
decision  must  be  made  as  to  whether  to  terminate  or  rephrase 
the  question  statement.   If  the  latter  choice  is  made,  the 
rephrased  question  will  be  included  In  the  following  rounds 
and  will  be  dropped  when  it  passes  the  test  of  stability  and 
when  no  further  benefits  can  be  expected  from  rephrasing  it 
(DeJanI  et  al .  .  1979). 
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For  example.  In  a  Delphi  study  with  20  respondents, 
consensus  would  be  achieved  with  unanimity  among  the  20, 
majority  with  11-19  responding  the  same  way,  bipolarlty  with 
a  10-10  split,  plurality  with  the  largest  subgroup  of 
respondents  between  2  and  9,  and  disagreement  with 
respondent  subgroups  numbering  0-1  (Dejani  et  al.,  1979). 
A  graphical  representation  of  Dejanl's  stopping  criterion  is 
presented  in  Figure  3. 

Upon  completion  of  each  Delphi  round,  the  mode,  mean, 
and  median  value  of  the  responses  can  be  calculated.  While 
It  would  be  desirable  to  obtain  100  %   consensus  from  all 
panel  members.  It  Is  extremely  unlikely.  Therefore,  it  Is 
necessary  to  determine  a  means  of  stopping  the  Delphi 
process  and  determining  a  response/value. 

For  the  purposes  of  this  study,  the  mode  value  for 
risk  was  used.   In  other  words,  a  cell  In  the  matrix  was 
assigned  a  risk  value  If  consensus  or  majority  of  the 
respondents  answered  the  same  way.   In  all  but  19  cells  (out 
of  156)  the  use  of  the  mode  was  an  acceptable  means  of 
assigning  risk  to  a  particular  activity  on  a  given 
slope/geology.  For  cases  where  using  the  mode  value  did  not 
result  In  one  value,  the  mean  of  all  responses  was  used  to 
assign  the  risk  value.  The  reader  must  bear  In  mind,  that 
the  matrix,  as  presented  In  this  thesis.  Is  not  as  yet 
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complete.  The  Delphi  procedure  Is  ongoing  and  these  results 
where  obtained  after  two  rounds.  Ceils  where  disagreement 
still  remained  after  two  rounds  are  delineated  by  an 
asterisk  (*).  The  final  erosion-Impact  matrix  after  2 
rounds  of  Delphi  Is  presented  In  Figure  4. 

DEFINITION  AND  JUST  I F ICAT iON(  FORf MATR IX-  FORMAT 

Design  of  a  simple  and  easy-to-use  erosion-Impact 
matrix  Is  essential.  The  matrix  developed  in  this  study  was 
patterned  after  the  one  used  in  Oregon  (Rlckert  et  ai., 
1978),  though  it  has  been  modified  to  better  fit  common 
management  concerns  in  western  Montana.  The  matrix  Is 
composed  of  an  X-axis  and  a  Y-axls.  These  axes  Include 
slope/geology  and  variable  land-use  activity 
classifications  (See  Figures  1  &  4). 

The  scope  of  this  research  deals  specifically  with  the 
effects  of  forest  practices  on  water  quality  degradation 
from  nonpoint  source  pollution.  These  activities  may 
Include  such  things  as  si  I v icui tura I  system  design,  logging 
methods,  site  preparation,  road  construction  and  road 
maintenance.  Therefore,  the  X-axis  of  the  matrix  is 
comprised  of  variable  land-use  activities  associated  with 
tree  harvest  operations. 
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The  X-axIs  Is  broken  down  Into  4  major  categories. 
These  are  as  follows:  1)  clearcuts;  2)  partial  cuts;  3)  site 
preparation;  and  4)  roads.  Each  category  Is  subsequently 
broken  into  different  practices  at  a  lower  level.  Clearcuts 
and  partial  cuts  are  subdivided  by  the  type  of  logging 
activities  which  are  possible  during  tree  removal.  Each 
category  delineates  harvest  operations  into  4 
sub-categories,  these  are:  1)  tracked  equipment,  2)  rubber 
tire  skldder,  3)  partial  suspension  of  logs  during  harvest, 
and  4}  full  suspension.   It  Is  necessary  to  distinguish 
among  logging  methods  due  to  the  variable  effects  of 
different  types  of  equipment  on  soil  stability, 
permeability,  and  Infiltration.  For  example,  a  log  that  Is 
skidded  to  a  landing  via  D-7  dozer  will  have  much  more 
Impact  on  soils  than  one  which  has  been  fully  suspended 
during  harvest  operations. 

Site  preparation  considers  3  basic  treatments: 
1)  Machine  piling/scarification;  2)  broadcast  burning;  and 
3)  other  methods  (hand  piling,  spot  burning,  chemical 
application  and  none).  Extensive  soil  damage  may  occur 
during  site  preparation.  Therefore,  It  is  necessary  to 
consider  the  potential  Impact  of  site  preparation  within  the 
matr  I  x . 
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Temporary  and  permanent  roads  are  considered  by  the 
process.   Megahan  and  Kldd  (1972),  studied  erosion  following 
road  construction  and  Jammer  logging  in  granitic  soils  of 
south-central  Idaho.   He  reported  an  Increase  In  erosion  due 
to  roads  of  770  times  that  of  natural,  and  220  times  natural 
due  to  logging  activity.  This  fact  Indicates  the  Importance 
of  roads  in  the  production  of  NPS  pollution  and  lends 
support  for  the  Inclusion  of  roads  into  the  erosion-Impact 
matr Ix. 

The  vertical  component,  or  Y-axIs,  of  the  Montana 
erosional-lmpact  matrix  consists  of  slope/soil  erodibiiity 
risk  classes.  The  first  part  of  the  Y-axIs  are  the  slope 
categories.   Initially,  four  slope  classes  were  defined 
based  upon  slope  breaks  recognized  as  limitations  with 
respect  to  land  management  activities  as  well  as  ease  In 
mapping.  These  are  outlined  as  follows: 

SLOPE  CLASSES  JUSTIFICATION 

0-20%  Used  in  the  Flathead  N.F. 

land  system  Inventory 
(Basko  et  al . ,  1983).   20% 
slope  is  recognized  as  the 
maximum  slope  for  which 
the  the  Universal  Soil 
Loss  Equation  may  be  used 
(Brown  III  et  al  .  ,1979) . 
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20-40%  40%  is  the  recommended 

upper  limit  for  certain 
mechanical  cultivation 
activities  (Basico  et  al. 
,1983). 

40-60%  60%  Is  recognized  as  the 

approximate  angle  of 
repose  for  sol  Is  and 
angular  rocic  fragments 
(Rlckert  et  a  I., 1978  and 
Basico  et  al  . ,  1983) 

>60% 

However,  after  the  1st  round  of  Delphi,  panel  members 
decided  that  the  first  slope  class,  0-20%,  was  too  broad  and 
as  a  result  of  committee  deliberation,  was  changed.   The  new 
slope  categories  for  the  erosion-impact  matrix  were  modified 
as  follows:  1)  0-5%:  2)  5-20%;  3)  20-40%;  and  4)  >40%. 
Justification  for  the  brealcs  at  20%  and  40%  are  the  same  as 
those  for  the  Initial  slope  classes. 

The  second  major  components  of  the  Y-axis  are  the  soil 
erodlbillty  classes.  There  are  3  proposed  soli  erosion 
classes.   They  are  1)  H-HIgh;  2)  M-Medlum;  and  3)  L-Low. 
These  classes  were  based  primarily  upon  groupings  of 
geologic  erosion  factors  for  soils  derived  from  various 
parent  materials  as  described  In  the  R1-R4  Sediment  Yield 
Prediction  Procedure  (Ciine  et  al.,1981).   The  categories 
are  reproduced  from  R1-R4  sediment  yield  models  In  Table  1, 

There  are  8  classes  of  materials  Identified  by  Andre 
and  Anderson  (1961)  described  In  Table  1.   These  8  classes 
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CATEGORIES 

GEOL.  EROSION  FACTOR 

L  Highly  Erodible 

a.  Add  Igneous  (Granitics) 

b.  Alliivium 

LO 

LOS 

2.  Moderately  Erodible 

a.  Schist 

b.  Soft  Sediments 
a  Hard  Sediments 

.75 
.66 
.52 

3.  Slightly  Erodible 

a.  Basiclgneous 

b.  Metamorpfaics  (Belt  Series) 

c.  Serpentine 

42 
.39 
.35 

TABLE  L  GEOLOGIC  EROSION  FACTORS  FOR  SOILS 
OF  VARIOUS  ROCK  TYPES  (Cline  et  aL,1981) 


are  based  upon  mean  surface  areas  (cm  )  of  silt  and  clay 
size  particles  over  the  mass  of  aggragated  silt  and  clay. 
Their  study  of  materials  In  northern  California  Indicated  a 
strong  relationship  between  these  values  and  erodlblllty 
(Cllne  et  al.,  1981).  Therefore,  since  they  are  in 
widespread  use  In  the  R1-R4  Sediment  Yield  model  and  have 
acquired  some  professional  acceptance.  It  was  decided  that 
use  In  the  Montana  eros lonai-impact  matrix  was  Justified. 
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CONCLUSION 

Risk  analysis  can  provide  resource  managers  with  a 
powerful  tool  to  assist  in  land-use  decision  making.  Recent 
laws  requiring  states  to  plan  for,  assess,  and  monitor 
waters  vulnerable  to  nonpoint  source  pollution,  have 
required  the  development  of  new  techniques  which  can  assist 
In  this  endeavor.  Modeling  watershed  responses  to  human 
encroachment  Is  undoubtedly  the  most  cost  effective  means  of 
predicting  possible  responses. 

This  study  was  designed  to  develop  a  methodology 
to  assist  In  assessing  NPS  pollution  and  cumulative 
watershed  effects.   This  was  done  by  combining  risk 
assessment  with  advanced  computer  technology  through  the  use 
of  geographic  Information  systems. 

Step  1  of  the  process  required  the  construction  of  a 
risk  matrix.   Associated  risk  ratings  were  assigned  to  the 
matrix  based  upon  proposed  si  I vicul tural  land-use  activity 
and  slope/geology  upon  which  It  was  to  take  place.  Values 
for  the  risk  matrix  were  obtained  through  a  variation  of  the 
Delphi  technique.  The  final  Montana  erosion-Impact  matrix 
for  risk  Is  presented  In  Figure  4. 

The  final  phase  In  this  procedure  Is  to  test  the  risk 
matrix  on  one  or  more  pilot  watersheds  In  western 
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Montana,  and  to  conduct  all  spatial  analysis  through  the  use 
of  GIS.  Presentation  of  GIS  risk  Integration,  and  results  of 
areal  calculations  for  the  2  pilot  studies  are  presented  In 
Chapter  2. 
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CHAPTER  2 
GIS  RISK  ASSESSMENT  APPLICATION 


INTRODUCTION 

Inadequate  processing  power.  Insufficient  memory 
capacity,  and  other  limitations  prevented  Implementation  of 
GIS  software  on  the  first  generation  of  microcomputers 
(Cooney  &  Tucker,  1986).   However,  In  recent  years  advances 
In  microcomputer  technology  have  made  the  use  of  GIS  by  land 
management  agencies  possible.  Because  environmental 
Information  Is  often  spatial,  geographic  Information  systems 
are  becoming  recognized  as  powerful  tools  for  resource 
management  (Robinson  et  al.  1987,  Bailey  1988,  Berry  and 
Sal  lor,  1987), 
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With  the  use  of  GIS,  planners  can  correlate  land  cover 
and  topographic  data  of  drainage  networks,  drainage  basin 
area,  and  terrain  configuration  (Walsh,  1985).  This  makes 
GIS  available  for  assessing  the  potential  effects  of 
land-use  activities  on  water  quality. 

Development  and  use  of  an  automated  GIS  can  expedite 
data  integration  problems  and  the  time-consuming  process  of 
synthesizing  tremendous  amounts  of  Information  for  spatial 
examination  of  nonpoint  pollution  (Walsh,  1985). 
Wilson  and  Thomas  (1977)  maintains  that  an  automated 
Information  system  through  which  geographically-referenced 
data  can  be  entered,  manipulated,  and  analyzed,  can 
Immeasurably  Improve  the  decision  making  process  of  land 
management  organizations. 

As  a  result,  many  federal  agencies  have  incorporated  or 
plan  to  Incorporate  GIS  Into  their  management  schemes.   For 
example,  the  Bureau  of  Land  Management,  which  has  management 
responsibility  for  over  340  million  acres.  Is  committed  to 
using  GIS  technology  to  help  manage  mineral,  range, 
wildlife,  and  forest  resources,  by  assessing  the  Impact  of 
proposed  development  (Hatch  1986,  Parker  1986). 

The  U.S.  Fish  and  Wildlife  Service  has  been  developing 
GIS  for  natural  resource  problems,  and  now  has  predictive 
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techniques  that  model  the  effects  of  channel  deepening  on 
the  distribution  of  various  shellfish  and  finflsh  species 
(Hatch  1986,  Robinson  et  al.  1987). 

Another  example  of  GIS  use  In  the  federal  government  Is 
the  U.S.  National  Park  Service.  The  Park  Service  created  a 
GIS  field  unit,  charged  with  developing,  applying,  and 
supporting  GIS  technology  throughout  the  Park  Service 
(Fleet,  1986).  A  prototype  GIS  has  recently  been  evaluated 
for  planning  aerial  spray  block  layout  for  fighting  spruce 
budworm  (Jordan  and  VIetlnghoff  1987,  Robinson  et  al., 
1987). 

ADVANTAGES/DISADVANTAGES  OF  GIS 


A  GIS  can  generally  be  described  as  a  computer  system 
for  entering,  storing,  managing,  retrieving,  transforming, 
analyzing,  and  displaying  spatial  data  (Cough  Ian  and  01  I  Iff. 
1988,  Cowan  1987,  Robinson  et  al.  1987).  However,  as  with 
any  advanced  technology  there  are  several  advantages  and 
disadvantages  that  must  be  considered  before  Implementation. 

The  major  strength  of  GIS  In  resource  management  Is  Its 
ability  to  link  spatial  and  nonspatlal  data  simultaneously. 
The  big  advantage  here  Is  that  large  amounts  of  nonspatlal 
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data  can  be  associated  with  graphical  representations 
(maps)  Immediately.   Information  such  as  stand  density, 
soil  type,  habitat  type,  erosion  classification,  ownership, 
and  land-use  history  can  be  coupled  to  spatial 
representations  such  as  polygons,  points,  or  lines. 
Nonspatlal  attributes  can  be  stored  and  accessed  through  a 
database  management  system  that  Is  either  externally  or 
Internally  I  inlced  to  the  GIS.   This  allows  for  Increased 
speed  and  efficiency  In  data  manipulation  and  may  expedite 
the  decision-making  process. 

Land-use  planners,  and  resource  specialists  use  maps  or 
Images  displaying  soil  type,  geology,  land-use,  hydrology, 
and  other  natural  resource  data  on  a  dally  basis.  Walsh 
(1985)  suggests  that  another  major  advantage  of  archiving 
these  data  in  a  GIS  Include  ease  of  retrieval,  variety  of 
output  products  to  fit  almost  any  need,  and  the  ability  to 
discover  and  display  Information  gained  by  testing  the 
interactions  between  natural  resources  phenomena  and  to 
organize  and  appraise  variable  coefficients  for  predictive 
mode  Is. 

Numerous  models  have  been  developed  to  predict 
hydrologic  watershed  responses  to  various  activities. 
Watershed-type  models  require  updates  of  land  cover, 
precipitation,  and  other  data  and  an  assessment  of  the 
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dynamic  spatial  patterns  affecting  runoff,  inf i I trat lon.and 
potential  nonpolnt  source  pollution  (Walsh,  1985).  The  use 
of  GIS  for  data  organization  and  manipulation  would  greatly 
benefit  model  speed  and  efficiency. 

Another  Intriguing  aspect  of  GIS  Is  the  potential  to 
receive  and  utilize  satellite  data.  Remote  sensing  is  a 
discipline  that  generates  volumes  of  spectral  data  on 
landscape  features  (Cough  I  an  and  Oil  Iff,  1988).  Various  GIS 
systems  are  able  to  Integrate  this  Information  for  use  in 
creating  quality  graphical  representations.  The  ability  of 
a  GIS  to  Integrate  raster  data  (grid  cell)  from  a  satellite 
sensor  with  data  digitized  from  a  map  sheet  Is  one  of  Its 
primary  strengths  (Robinson  et  al.,  1987). 

On  the  other  hand,  there  are  several  disadvantages  to 
be  considered  before  Implementing  a  GIS.  The  first  major 
consideration  is  money,   implementing  a  GIS  requires  a 
substantial  Initial  Investment,  both  in  dollars  and  In 
time  to  digitize  maps,  build  a  database  and  develop 
customized  reporting  (Devlne  &  Field  1986,  Relslnger  &  Davis 
1987).   Software  modules  alone  range  from  $10,000-$100,000 
for  micro-based  systems  (Cough  Ian  and  01  I  iff,  1988),  and  the 
cost  of  keeping  data  current  can  quickly  outstrip  the 
Initial  Investment  (Devlne  &  Field  1986). 
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AntenuccI  (1986)  recommends  a  detailed  feasibility 
cost-benefit  analysis  be  conducted  before  Implementation. 
He  maintains  that  the  feasibility  analysis  will  allow  an 
organization  to  determine  If  conversion  to  an  automated  GIS 
can  be  Justified  by  calculating  the  costs  to  be  Incurred  by 
the  organization  compared  against  the  benefits  of 
automat  Ion. 

Four  additional  disadvantages  of  GIS  are:   1)  a 
substantial  amount  of  hardware  must  be  available.  Including 
a  digitizer  and  computer  plotter,  a  photogrammetr Ic 
tr langulat Ion  Instrument,  and  relatively  large  computer 
memory  capacity;  2)  although  automated  mapping  processes 
requires  less  total  time  than  traditional  methods,  longer 
lead  times  and  careful  planning  are  required  to  iinic  many 
Intermediate  taslcs  Into  a  completed  data  base;  3)  GIS  do 
not  lend  themselves  to  single-purpose  operations.   In  other 
words,  they  are  most  efficient  when  creating  a  mult  I  faceted 
information  base  where  the  Integration  of  several  separate 
data  sources  Is  required.  And  finally,  4)  it  Is  necessary 
that  people  from  different  disciplines  woric  together  In 
building  and  using  a  common  database  (Martin,  1985). 


Page  38 

The  final  disadvantage  can  be  the  advanced  and 
Intricate  nature  of  CIS  technology  Itself.  Most  software 
and  hardware  packages  are  complicated  and  can  be  very 
difficult  to  understand.  Therefore,  staffing  Is  a  primary 
consideration.  Adequately  trained  personnel  are  a  must  If 
maximum  benefit  and  equipment  potential  Is  to  be  realized. 
However,  GIS  technology  can  be  very  difficult  to  teach 
(Berry,  1986).  Practical  experience  Is  required  as  well  as 
theory,  yet  very  few  classrooms  can  provide  extensive 
hands-on  learning.  More  systems  have  failed  due  to 
inadequately  trained  staff,  than  for  any  other  reason 
(AntenuccI ,  1986) . 

PREVIOUS  USE  OF  GIS  IN' WATERSHED  MANAGEMENT 

While  there  are  some  negative  Implications  that  may 
arise  In  the  Implementation  of  a  new  GIS,  the  benefits  of 
automated  mapping  can  be  Immeasurable.   There  are  numerous 
examples  where  GIS  have  been  used  for  natural  resource 
management  planning,  specifically  for  watershed  management 
and  NPS  pollution  control. 

Atkinson  (1987)  conducted  a  study  on  21  watersheds 
draining  Into  the  Trinity  River  in  the  Dallas-Fort  Worth 
(Texas)  Metroplex.  The  study  utilized  geographic 
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Information  systems  to  address  the  question  of  whether  or 
not  Trinity  River  water  quality  could  be  protected  from  NPS 
pollution  by  using  wet  detention  basins.  Results  Indicated 
that  a  GIS  modeling  approach  was  successful  In  determining 
the  most  feasible  sites  for  the  detention  basins  and  was 
Instrumental  In  the  planning  and  design  process. 

A  similar  study  by  Berry  and  Sailor  (1987)  describe  the 
use  of  generalized  GIS  for  analyzing  the  spatial  aspects  of 
storm  runoff  prediction  using  the  US  Soil  Conservation 
Service  (SCS)  technique  (Sailor  &  Berry,  1980).  The 
objective  of  their  study  was  to  familiarize  environmental 
planners  with  application  of  geographic  Information  system 
to  storm  runoff  prediction.   The  results  of  their  study 
Indicated  that  the  use  of  GIS  to  generate  watershed  data  to 
be  used  In  the  SCS  storm  runoff  prediction  method  greatly 
enhanced  precision  and  speed  over  manual  procedures.  The 
management  and  analytical  capabilities  of  GIS  technology 
provided  them  with  the  opportunity  to  fully  Integrate 
spatial  conditions  Into  hydrologic  inquiries.   They 
concluded  that  with  advent  of  GIS  systems  for  microcomputer 
environments  and  the  Increasing  availability  of 
digitally-mapped  data,  GIS  technology  will  play  an  important 
role  In  other  spatial  studies  In  environmental  planning,  as 
well  as  in  storm  water  monitoring  (Sailor  and  Berry,  1980). 
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Prato  et  al.  (1989),  described  how  geographic 
Information  systems  were  used  to  assemble  and  retrieve  the 
physical  parameters  required  to  estimate  sheet  and  rill 
erosion  and  water  quality  Impacts  of  multiple  resource 
management  techniques  on  16  farms  in  northern  Idaho.   One  of 
the  study's^maln  objectives  was  to  demonstrate  how  a  GIS 
could  be  used  to  estimate  soil  erosion  from  fields  and  farms 
In  a  watershed  and  to  estimate  sediment  and  nutrient 
pollution  In  receiving  waters  for  alternative  resource 
management  (Prato  et  al.,  1989).   Results  indicated  that  the 
use  of  GIS  in  data  acquisition  and  spatial  analysis  was 
paramount  In  the  overall  success  of  the  study. 

OBJECTIVES  (GIS) 

As  a  result  of  recent  successes  In  the  application  of 
GIS  in  natural  resource  management,  particularly  in 
watershed  management,  the  integration  of  GIS  In  the 
assimilation  and  development  of  a  risk  analysis  procedure 
for  use  In  forested  lands  is  presented.   This  procedure  was 
tested  on  2  pilot  watersheds  In  western  Montana.   The  first 
pilot  study  was  conducted  on  Howard  Creek,  a  third  order 
watershed  on  the  Lolo  National  Forest.  The  second  pilot 
watershed,  Jones  Meadow  Creek,  is  in  the  Lubrecht 
Experimental  Forest. 
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The  Study  Areas 

Howard  Creek  is  located  approximately  32  kilometers 
southwest  of  Missoula,  Montana,  In  the  Lolo  National  Forest. 
It  was  selected  as  a  pilot  study  area  because  of  its 
diversity  in  past  land-use  history.  Currently,  Howard  Creek 
Is  managed  by  the  Lolo  National  Forest,  Champion 
International,  and  Plum  Creek  Timber,  Inc. 

Jones  Meadow  Creek  Is  located  in  the  Lubrecht 
Experimental  Forest,  managed  by  the  University  of  Montana, 
School  of  Forestry,  It  was  selected  as  a  pilot  study  area 
because  most  of  the  spatial  Information  required  to 
calculate  risk  values  had  been  previously  digitized. 

Topography 

The  Howard  Creek  watershed  encompasses  an  area  of 
approximately  5015  hectares.   It  exhibits  a  dendritic 
drainage  pattern  typical  of  small  watersheds  In  western 
Montana,   Relief  within  the  basin  ranges  from  a  maximum 
elevation  of  1767.8  m  to  a  low  of  1188,7  m;  with  an  average 
watershed  elevation  of  approximately  1480  m.  The  main  stem 
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of  Howard  Creek  Is  oriented  east-west  while  the  3 
sub-basins,  Teppee ,  North  Fork,  and  Krystal,  are  oriented 
north-south  (See  Figure  5). 

The  Jones  Meadow  watershed  totals  635.4  hectares. 
Relief  within  the  basin  ranges  from  1250  m  to  1615  m.   The 
watershed  Is  oriented  north-south  with  primary  aspects 
facing  north. 

CI Imate 

The  Lolo  National  Forest  Is  under  a  modified  maritime 
climate  regime  (Saslch  and  Lamotte-Hagen,  1989).   The 
Continental  Divide  creates  a  physical  barrier  which  effects 
climatic  processes  In  Montana.  Areas  west  of  the  Divide  are 
dominated  by  a  maritime  (North  Pacific  Coast)  climate  and 
those  east  of  the  Divide  are  dominated  by  a  continental 
climate.  Howard  Creek  Is  located  west  of  the  Continental 
D I V I de . 

Precipitation  results  from  orographic  and  frontal 
activity  associated  with  low  pressure  systems  originating 
off  the  Pacific  coast.  Approximately  two-thirds  of  the 
precipitation  received  falls  as  snow.   Howard  Creek  averages 
100  cm  of  annual  precipitation.   Temperatures  range  from  an 
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average  of  -5.5  C  In  January  to  an  average  of  19.4  C  In 
July.   Estimates  for  Howard  Creek  can  be  extrapolated  from 
Information  compiled  at  the  nearest  weather  station  located 
at  Lolo  Hot  Springs.   A  graph  of  the  average  annual 
precipitation  for  Lolo  Hot  Springs  Is  presented  in  Figure  6. 

Climate  for  the  Lubrecht  Experimental  Forest  Is 
described  as  modified  temperate  continental  regime. 
Modified  temperatures  result  from  maritime  influences 
originating  In  the  northern  Pacific  ocean.  Monthly  average 
temperatures  for  Lubrecht  Experimental  Forest  headquarters, 
approximately  1.5  kilometers  (km)  to  the  west,  can  be  used 
to  estimate  climate  for  the  Jones  Meadow  watershed. 
Temperature  extremes  at  the  headquarters  facility  ranged 
from  a  summer  high  of  40  C  to  a  summer  low  of  -6.0  C; 
winter  extremes  range  from  a  minimum  of  -42.5  C  to  a 
maximum  of  20  C  (Goetz,  1989). 

Average  annual  precipitation  for  the  Lubrecht 
experimental  watershed  Is  about  45.5  cm.   Approximately  44% 
falls  as  snow  in  the  winter  months  (Nov-Mar)  and  24%  fails 
during  the  summer  (Jun-Aug). 
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Figure  6:  Annual  Precipitation  and  Snow  Data 
For  Lolo  Hot  Springs  (1959-1984) 
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Geology 

Soil  parent  materials  are  In  part  determined  by 
bedrock  and  surflcial  geology.  The  Lolo  National  Forest 
Land  Systems  Inventory  (Saslch  and  Lamotte-Hagen,  1989) 
maps  the  primary  geologic  groups  found  In  Howard  Creek  as 
metased imentary  and  undifferentiated. 

Metasedlmentary  rocks  contain  parent  materials  derived 
from  Belt  Super  Group  quartzite,  argil  lite,  and  siltlte. 
Rock  fragment  hardness  is  variable  depending  upon  the  degree 
of  rock  weathering.  Weathering  Is  dependent  on  associated 
faults,  preponderance  of  arglllltes,  and  calcium  carbonate 
content.  These  materials  were  classified  as  either  L-low  or 
M-moderate  within  of  the  erosion-Impact  matrix. 

Undifferentiated  geology  is  composed  of  materials 
derived  from  Belt  Super  Group  metasedlmentary  rocks  or 
weakly  weathered  granitic  rocks.  Materials  Include  alluvium 
on  terraces  and  flood  plains;  shallow  soils  on  flood  scoured 
foot  slopes  and  stream  breaklands,  strongly  frost  churned 
broadly  convex  ridges,  and  glacial  outwash  on  plains  (Sasich 
and  Lamotte-Hagen,  1989).  This  group  would  be  classified  as 
H-hlghly  erodible  in  the  Montana  erosion-Impact  matrix. 


Page  47 

Geology  for  Howard  Creek  was  mapped  using  GIS.  A 
graphic  representation  of  the  3  erosion  hazard  groups,  high, 
medium,  and  low,  can  be  seen  In  Figure  7. 

Brenner  (1964)  mapped  geology  of  Lubrecht  and  this  Is 
summarized  by  NImlos  (1986)  In,  Sol ls0of0Lubrecht 
Exper lmentai0Forest  .  These  include  belt  rocks,  limestone, 
granite,  tertiary,  and  unconsolidated  transported  material. 

The  first  type  of  parent  material  found  in  Jones  Meadow 
are  belt  rocks.  They  are  the  oldest  rocks  In  the  Beit  Super 
Group  and  were  deposited  during  the  precambrian  about  one 
billion  years  ago.  They  are  originally  formed  as  marine 
deposits  which  have  metamorphosed  into  quartzlte,  argll.lite, 
and  siltites  (NImlos,  1986).  These  deposits  are  included  in 
the  erosion-Impact  matrix  In  the  LOW  erodlbility 
classification. 

The  next  major  parent  material  found  In  the  Jones 
Meadow  watershed  are  tertiary  deposits.  Tertiary  rocks  are 
derived  from  sediments  deposited  60  million  years  ago. 
These  became  consolidated  in  weakly  cemented  and 
interbedded  conglomerate,  sandstone,  siltstone,  and  mudstone 
(NImlos,  1986).  Tertiary  parent  materials  are  considered 
with  alluvium  as  highly  erodible  and  are  assigned  a  HIGH 
classification  in  the  erosion-impact  matrix.  Figure  8  shows 
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Erodibilitv  Classification 


Figure  8:   Geology  Map  for  Jones   Meadow 
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the  distribution  of  erosion  classification  groups  for  Jones 
Meadow. 

Land-Use  History 

Ownership  of  Howard  Creek  Is  divided  among  three 
entitles:  1)  the  Lolo  National  Forest;  2)  Champion 
International;  and  3)  Plum  Creek  Timber  Inc.   Ownership  was 
mapped  using  GIS  and  can  be  seen  In  Figure  9.   Champion 
International  controls  approximately  34%,  the  USFS  owns 
approximately  42%,  and  the  remaining  23%  Is  owned  and 
managed  by  Plum  Creek  Timber  Inc. 

Timber  harvest  activities  In  recent  years  have  Impacted 
27.4%  of  Howard  Creek's  total  acreage.   In  a  recent  report 
to  the  forest  supervisor,  cumulative  watershed  effects  have 
been  Judged  as  extremely  detrimental  (Munther  et  al.,  1987). 
In  fact,  land-use  activities  In  Howard  Creek  have  resulted 
In  a  current  sediment  load  Increase  of  50%  and  current  water 
yield  Increase  of  8%. 

Lubrecht  Experimental  Forest  Is  on  state  owned  lands  and 
Is  managed  by  the  School  of  Forestry,  University  of  Montana. 
It  consists  of  a  tract  of  11,331.5  hectares  and  is  located 
about  48  kilometers  (km)  northeast  of  Missoula.   The  forest 
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was  established  In  1937  to  foster  research  on  natural 
resources,  but  It  Is  also  used  for  forestry  education  and 
demonstration  of  forest  management  techniques  (NImlos, 
1986). 

Land-use  activities  In  the  Jones  Meadow  watershed  have 
not.  In  terms  of  silviculture,  been  extensive.   Only  seven 
harvest  blocks  have  been  cut  In  the  past  decade  and  these 
were  small  by  comparison.  The  area  was  subjected  to 
experimental  treatments  that  had  minimal  Impact  on  the 
stability  of  soil  and  other  natural  resources. 

GIS  Processes 

In  1988  the  University  of  Montana,  School  of  Forestry, 
began  the  process  of  acquiring  a  Geographic  Information 
System  (GIS).  Choosing  a  GIS  package  that  was  cost 
efficient  and  whose  foundation  focused  on  natural  resource 
management  were  essential.  The  School  of  Forestry  selected 

PAKCAP. 

PAMAP  was  developed  by  PAMAP  Graphics  LTD  In  Victoria, 
British  Columbia.   it  is  one  of  the  most  powerful  analytic 
tools  available  in  today's  automated  mapping  environment 
(PAMAP.  1989). 

The  need  for  computer  mapping  and  data  manipulation 
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spawned  from  the  ever-growing  value  of  information  generated 
by  such  techniques  as  remote  sensing.  Paper  maps,  simply 
could  not  keep  up  and  are  extremely  limited  In  what  can  be 
displayed  visually.   In  addition  to  performing  as  poor 
databases,  paper  maps  were  also  Ill-suited  to  depict  the 
results  of  various  analytical  functions. 

As  a  result,  PAMAP  was  Instrumental  In  creating  the 
ability  to  merge  computers  and  cartography.   PAMAP  was 
developed  In  the  early  I980's  and  soon  proved  Itself  as  a 
mapping  tool  ideally  suited  to  processing  and  analyzing 
great  volumes  of  data,  and  producing  high-quality  maps 
(PAMAP,  1989). 

Its  most  valuable  asset  is  its  ability  to  model  the 
physical  world.   Since  data  can  be  accessed,  transformed, 
and  manipulated  Interactively,  GIS  allows  one  to  test 
processes,  analyze  trends,  and  anticipate  the  outcome  of 
planning  Initiatives  without  ever  touching  whatever  part  of 
the  physical  world  being  considered  (PAMAP,  1989). 

PAMAP  software  consists  of  4  primary  modules  which 
perform  5  basic  functions.  The  modules  include:  1)  MAPPER; 
2)  TOPOGRAPHER;  3)  ANALYZER;  and  4)  INTERPRETER.   The  5 
fundamental  functions  are:  1)  Data  input  and  Verification  - 
which  covers  all  aspects  of  transforming  data  from  existing 
maps,  field  measurements,  and  satellite  sensors  into  a 
usable  digital  format;  2)  Data  Storage  and  Database 
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Management  -  defines  the  structure  and  organization  of  the 
data  map  position.  Its  corrections  to  other  elements 
(topology)  and  the  non-graphic  attributes  of  geographical 
elements  such  as  points,  lines,  and  area;  3)  Data  Output  and 
Presentation  -  output  is  directed  to  either  a  plotter  or 
color  copier;  4)  Data  Transformation  Functions  -  Including 
error  removal  and  data  analysis  that  can  affect  the  location 
of  the  data.  Its  non-spatial  attributes,  or  both  ;  and  5) 
Interaction  With  User  -  allows  direct  access  to  systems  data 
and  lets  the  user  manipulate  attributes  If  necessary  (PAMAP, 
1989). 

MAPPER 

GIS  MAPPER  Is  PAMAP's  Interactive  data  entry  module. 
It  Is  used  to  digitize  new  maps  and  edit  existing  ones. 
MAPPER  requires  that  paper  maps  be  converted  Into  digital 
form.  This  can  be  done  manually,  or  with  a  optical 
digitizing  scanner.  Once  the  map  data  Is  In  the  central 
processing  unit  (CPU),  the  full  range  of  GIS  tools  Is 
accessible  by  the  user.  MAPPER  can  easily  modify  existing 
lines,  and  redraw  them  if  needed. 

MAPPER  can  Input  data  in  four  coordinate  systems:  1) 
Geographic  Units  (Longitude  and  Latitude),  2)  UTM 
projection,  3)  polyconic  projection,  and  4)  Lambert 
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conformed  projection.  Additionally,  each  map  can  Include  64 
layers,  or  levels,  of  geographic  Information,  each  supported 
by  attribute  data. 

ANALYZER 

ANALYZER  opens  the  door  to  PAMAP's  powerful  analysis 
tools  by  converting  the  lines,  points  and  polygons  of  map 
line  work,  together  with  attribute  data,  Into  a  grid-cell 
format.  As  it  converts  a  thematic  map  layer  to  grid-ceil 
format,  ANALYZER  automatically  forms  polygons,  computing 
their  area  and  perimeters  and  storing  them  In  the  data  base. 

ANALYZER  analyses  geographic  data  by  performing  2 
taslcs:  conversion  and  overlaying.  ANALYZER  may  simply 
convert  one  map  level  from  a  given  format  to  another,  and 
then  use  the  resulting  map  for  overlaying  operations. 
Overlaying,  requires  one  raster  cover  to  be  overlayed  upon 
another  raster  cover.  This  procedure  creates  a  third  raster 
map  for  which  resulting  statistics  can  be  generated. 

TOPOGRAPHER 


TOPOGRAPHER  Is  a  sophisticated,  yet  easy-to-use  tool 
for  manipulating  three-dimensional  data.   It  can  accept  and 
integrate  both  vector  (points)  or  raster  (grid)  formats. 
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TOPOGRAPHER  can  generate  elevation  models  based  on  a  variety 
of  Interpolation  routines.  Derivative  surfaces  such  as 
slope  and  aspect  may  be  produced,  as  well  as  perspective 
views  (PAMAP,  1989).  Another  very  Important  feature  of  this 
module  Is  its  ability  to  Import  elevation  and  data  from 
numerous  sources,  including  irregular  ASCII  data  from 
external  sources  and  contours  digitized  in  MAPPER. 

INTERPRETER 

The  final  PAMAP  module,  and  one  of  the  most  important, 
is  INTERPRETER.  PAMAP  Graphics  LTD  has  worlced  for  several 
years  with  the  Canada  Centre  for  Remote  Sensing;  recently 
developing  a  technique  for  integrating  the  two  technologies 
(PAMAP,  1989);  INTERPRETER  is  the  result. 

It's  major  function  Is  to  convert  remote  sensing  data 
Into  a  high-resolution  graphic  display  In  MAPPER.  This 
malces  INTERPRETER  one  of  the  most  valuable  features  of 
PAMAP's  GIS.   It  provides  a  means  for  resource  managers  to 
automatically  incorporate  the  hordes  of  data  generated  by 
satellite  platforms  Into  a  useable  format. 

All  modules  in  PAMAP  are  linked  to  one  another  through 
GIS  (Figure  12). 
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Data  Requirements 

In  order  to  apply  risk  analysis  techniques  on  any  given 
watershed.  Initial  data  assimilation  Is  required. 
Fundamentally,  Information  on  4  primary  spatial  attributes 
must  be  obtained  before  risk  analysis  and  digitizing  can 
begin.   These  are:  1)  sol  I /geo logy ;  2)  topography 
(elevation);  3)  land-use  history;  and  4)  system  and 
temporary  road  networks. 

Information  on  soils  can  be  obtained  from  geology  and 
soil  maps  that  are  produced  by  the  U.S  Geological  Survey, 
the  Soil  Conservation  Service,  U.S.  Forest  Service  Service 
and  from  various  state  agencies.   It  Is  Important  to  obtain 
accurate  soli  data  and  to  remain  consistent  In  Its 
Integration  into  the  GIS.  jytost  soil  maps  Include  soil 
classifications  based  on  soli  taxonomy.  Some  maps  will 
Include  information  about  parent  materials  from  which  soils 
are  derived.   For  application  of  the  Montana  erosion-impact 
matrix.  It  was  decided  that  for  ease  of  mapping  and  since 
the  R1/R4  Sediment  yield  model  Is  based  upon  them,  that 
parent  materials  and  not  soli  classifications  be  mapped. 

SCS  soli  maps  for  Howard  Creek  and  Jones  Meadow  were 
used  and  soil  classifications  were  grouped  into  parent 
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materials  categories  for  mapping  purposes.  The  resuits  are 
shown  as  geologic  (erodlbl I  I ty)  maps  In  Figure  7  and  8. 

The  next  vital  link  In  the  digitizing  process  Is 
topographic  Information  (contours).  Topographic  data  can 
pose  a  major  problem  in  terms  of  man  hours  and  complexity  In 
digitizing.   Fortunately,  topographic  information  for  most 
areas  Is  accessible  from  the  USGS  in  the  form  of  Digital 
Line  Graphs  (DLG)  or  a  Digital  Elevation  Model  (DEM).   This 
saves  numerous  hours  of  manual  digitizing  time.   However, 
expense  and  delays  in  shipping  may  make  DLGs  and  DEMs  not 
practical.   Therefore,  as  In  the  case  of  this  study,  manual 
digitizing  Is  often  the  answer.   It  took  approximately  24 
hours  at  the  digitizing  pad  to  input  contour  Information  for 
Howard  Creek. 

The  third  element  which  is  essential  In  the  analysis 
process  Is  past  land-use  history.  This  can  pose  a  major 
hurdle  In  terms  of  accuracy  and  timeliness  of  assimilation. 
Land-use  history  may  be  obtained  from  local  or  regional 
planning  offices,  federal  agencies  such  as  the  Forest 
Service,  Bureau  of  Land  Management,  or  the  USGS  and  from 
private  Industry  which  may  control  or  own  a  portion  of  the 
watershed.   If  the  basin  lies  in  a  relatively  pristine 
location,  as  does  Howard  Creek,  then  urbanization  is  not  a 
factor  and  acquiring  data  may  be  much  easier.   For  this 
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project,  all  Information  on  land-use  was  obtained  from  the 
primary  landowners.   Sources  of  data  ranged  from  paper  maps 
In  old  filing  cabinets  to  advanced  computer  graphics  and 
files.   Land-use  history  for  the  past  decade  was  obtained 
and  maps  showing  all  major  silviculture  operations  are  shown 
In  Figure  10  and  Figure  11,  Howard  Creek  and  Jones  Meadow 
respect Iveiy. 

The  final  geographic  spatial  attribute  needed  Is  the 
road  network.  Since  roads  can  have  major  negative 
effects  on  runoff  characteristics,  roads  are  vital  for 
inclusion  Into  the  matrix  design  and  In  risk  assessment 
(Megahan  1980,  Hasklns  1986,  and  Rickert  et  al.  1978). 

As  with  many  other  forms  of  spatial  Information,  road 
systems  can  be  found  on  paper  maps.   The  USGS  1:24000  scale 
maps  are  usually  outdated.  Therefore,  road  system 
information  must  be  obtained  from  the  various  landowners. 
Perhaps  the  very  best  source  for  acquiring  accurate  road 
information  is  through  aerial  photography.  This  was  done 
for  Howard  Creek  and  proved  to  be  Invaluable. 

Once  all  required  information  has  been  assimilated  It 
Is  necessary  to  consolidate  it  onto  one  or  two  databases. 
While  paper  maps  have  their  utility,  It  Is  not  a  good  Idea 
to  rely  upon  them  as  a  database  from  which  to  begin 
digitizing.   Paper  map  products  are  easily  bent,  torn. 
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stretched,  and  destroyed.   It  Is  strongly  suggested  that 
mylar  map  sheets  be  used  for  all  digitizing  purposes.  Mylar 
Is  a  plastic  material  made  from  a  synthetic  polyester 
compound,  polyethylene,  and  can  be  obtained  with  all  map 
information  printed  on  It  from  local  printing  or  blueprint 
shops.   However,  It  Is  very  expensive  in  comparison  to  paper 
map  products  but  will  be  well  worth  It  In  the  long  run. 

Once  ail  spatial  Information  Is  combined  onto  one  mylar 
map  sheet  digitizing  may  begin.   Hence,  starts  the  long  and 
somewhat  tedious  process  of  data  entry. 

Data  Input  To  GIS 

Data  entry  Into  any  GIS  Is  always  an  extremely 
challenging  task.  The  vast  amounts  of  spatial  Information 
for  even  the  most  basic  map  can  seem  Insurmountable. 
Patience,  diligence,  and  dedication  are  required  if  a 
project  Is  to  be  successful.  There  are  7  primary  steps 
involved  In  the  initial  data  Input  and  risk  assessment. 
These  steps  are  not  software  dependent.   They  Include:  1) 
GIS  preparation;  2)  level  definition;  3)  conversion  (vector 
to  raster);  4)  addition  of  non-spatial  data;  5)  data  file 
manipulation;  6)  three-dimensional  analysis;  and  7)  spatial 
analysis. 
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The  first  task  Is  GIS  preparation.   It  Is  not 
practical  to  begin  digitizing  without  extensive  forethought. 
Most  GIS  can  conform  to  user  specific  needs.   In  other 
words,  data  files  can  be  set  up  with  specific  requirements 
In  mind.   Definition  of  map  parameters  are  required  before 
data  can  be  entered. 

These  Include  such  things  as  map  scale,  UOR  (Units  of 
Resolution),  and  coordinate  systems.  For  example,  PAMAP  was 
designed  In  Canada  and  It  was  originally  developed  using  the 
metric  system  (Systems  International),  therefore  Howard 
Creek  was  configured  entirely  In  metric  units.  The  basic 
UOR  is  the  centimeter. 

, Another  Important  step  In  the  design  of  a  fluid  working 
system  Is  level  definition.   It  Is  essential  that  planning 
and  organization  In  terms  of  data  flies,  levels,  be 
maintained  from  the  onset.  There  are  64  levels  available  in 
PAMAP.   Each  level  can  Include  3  types  of  file:  1)  vector 
(points);  2)  polygonal  (grid-ceil);  and  3)  surface  (also 
grid-celi).  Therefore,  the  total  number  of  possible  levels, 
files.  Increases  to  192.   It  Is  Imperative  that  a  good 
Index,  or  reference  system,  be  used  to  minimize  confusion. 
This  problem  is  not  limited  solely  to  PAMAP  software,  but 
Is  un I  versa  I  . 
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Once  vector  data  has  been  put  In  the  GIS.  the  next 
major  step  Is  conversion.  The  power  of  spatial  analysis 
comes  from  the  ability  for  the  GIS  to  quickly  assimilate 
Information  areaily.  GIS  cannot  calculate  area  based  upon 
vector  data.   In  the  vector  to  raster  conversion  process, 
polygons  are  created  and  a  grid-cell  network  Is  produced 
based  upon  where  and  how  the  vector  Information  was  Input. 
This  takes  the  form  of  polygonal  or  surface  cover.  At  this 
point,  non-spat  lai  attribute  information  can  be  attached  to 
the  graphics  files  and  various  calculations  such  as  area  are 
possible.  The  raster  file  Itself  consists  of  a  series  of 
pixels  (cells)  that  represent  the  surface  of  the  ground  as 
depicted  on  the  map.  Pixel  size  can  be  adjusted  based  upon 
the  users  desired  resolution  and  Is  a  key  factor  In  the 
planning  process. 

The  next  major  step  In  developing  a  useable  product  can 
be  considered  the  most  Important.   It  Is  the  creation  of 
polygon  database  user  attributes,  or  more  precisely,  the 
link  of  non-spatial  data  to  digitized  spatial  graphic  files. 
Here  Is  where  the  GIS  gets  Its  tremendous  advantage  over 
traditional  cartographic  processes.  The  ability  to  store 
vast  amounts  of  non-spatial  Information  which  Is  immediately 
accessible,  makes  GIS  an  Invaluable  planning  tool.   Thought 
and  foresight  must  go  into  creating  databases.   Although  GIS 
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is  capable  of  storing  Incredible  amounts  of  Information,  It 
Is  best  to  keep  non-spat lal  Information  as  precise  as 
possible.  The  use  of  abbreviations  Is  suggested.  Once 
attribute  Information  is  entered.  It  is  possible  to  use 
"theming"  or  coloring  to  highlight  areas  on  the  map  that  are 
of  particular  Interest. 

Once  all  data  files  have  been  created  and  are  In  the 
GIS,  the  power  of  the  system  can  be  realized.  Levels  may  be 
overlaid,,  or  combined,  creating  a  resulting  map.   Files  can 
be  manipulated  to  generate  Information  with  optimal  land-use 
In  mind.   Corridor  calculations  can  lay  out  streams  I de 
management  zones  and  determine  acreage  of  lost  timber 
volume,  by  species  if  desired.  Area,  distance,  and 
perimeter  information  can  be  determined  Instantly. 
Polygonal  covers  can  be  combined  and  the  resulting  level  can 
be  used  for  analysis  and  planning. 

For  example,  2  polygonal  covers,  one  containing  species 
Information  and  one  containing  parent  materials,  are 
combined.   The  resulting  polygonal  cover  contains  both 
attributes.  The  user  can  query  the  system  and  obtain 
instant  results  on  the  new  cover.   If  Information  Is  desired 
on  all  lodgepole  pine  that  are  located  on  soils  derived  from 
colluvium,  the  GIS  will  instantly  highlight  all  areas 
meeting  the  specified  criteria. 
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The  sixth  step  Is  topographic  analysis.  Once  elevation 
Is  put  Into  a  data  file,  the  possibilities  for 
three-dimensional  manipulation  are  limitless.   Slope  and 
aspect  polygons  can  be  generated.  Perspective  views,  or 
viewablllty  from  a  certain  point  In  the  watershed.  Is  also 
possible.  This  factor  can  be  very  Important  for  planning 
purposes  where  visual  aesthetics  are  essential,  such  as  In 
timber  harvest  operations.  Furthermore,  detailed  analysis 
of  basin  relief  Is  possible.  This  Is  of  particular  Interest 
to  forest  hydrologlsts.   In  PAMAP,  all  3-D  analysis  Is 
generated  from  a  digital  elevation  model  (DEM).  A  DEM  for 
Howard  Creek  Is  presented  In  Figure  13. 

The  final  step  In  the  risic  analysis  process  and  In  data 
entry  Is  spatial  analysis.  Spatial  analysis  is  basically  an 
in  depth  examination  of  the  physical  geographic  features 
that  are  present  in  a  given  watershed.  Files  and  software 
can  be  manipulated  in  order  to  obtain  relevant  information. 
To  calculate  risic  it  is  necessary  to  determine  the 
percentage  of  the  watershed  that  is  sensitive  to  nonpoint 
source  pollution  production  or  further  land-use.  These 
areas  can  be  retrieved  from  the  database  and  risk  on  an 
areal  basis  can  be  calculated. 
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GIS  OUTPUT 

A  Geographic  Information  System  Is  a  computer  assisted 
map  making  technology.  The  primary  product  of  any  GIS  is  a 
graphic  representation.   In  other  words,  a  GIS  makes  maps. 
However,  information  such  as  area  and  distance  can  be 
calculated  and  retrieved  visually.  If  not  on  paper. 

For  the  purposes  of  this  study,  the  major  output  was 
In  the  form  of  areal  calculations  and  maps  of  high  risk 
areas.  The  area  of  each  watershed  classified  as  either 
H-HIgh,  M-Moderate,  L-Low,  In  terms  of  erodlblllty,  was 
determined  with  respect  to  geology  and  slope.  Previous 
land-use  activities  were  also  mapped  and  area  calculations 
were  generated. 

The  link  between  the  Montana  erosion-impact  matrix  and 
the  GIS  for  both  watersheds  in  terms  of  areas  sensitive  to 
NPS  pollution  production  Is  presented  In  the  following 
calculat Ions. 

CALCULATIONS 

There  are  four  basic  steps  Involved  In  total  watershed 
risk  calculations.  They  are:  1)  calculation  of  erosion 
classifications  by  area;  2)  calculation  of  risk  areal ly 
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based  on  geology  and  slope-,  3)  calculation  of  risk  areally 
based  on  past  and  planned  land-use  activities;  and  4) 
summation  of  risk  values  for  slope/geology  and  land-use 
act Ivl t les. 

Step  one  requires  that  areal  calculations  based  upon 
erosion  classifications  be  determined.  This  Is  relatively 
simple  In  that  the  GIS  can  Instantly  generate  the  total  area 
of  the  watershed  which  have  been  classified  as  H-HIgh 
erodibility.  M-Moderate  erodlbillty,  and  L-Low  erodlbiiity. 
At  this  point,  the  values  obtained  for  H,M,  and  L  land  areas 
are  divided  by  the  total  acreage  of  the  watershed.  This  was 
done  for  both  pilot  watersheds  and  is  presented  In  Table  3. 

Calculations  for  risk  based  upon  slope  and  geology  are 
presented  in  Table  11.  They  are  based  upon  areal 
distribution  of  erosion  classifications  as  determined  In 
step  1.  The  Montana  Erosion-Impact  Matrix  was  developed  on 
a  relative  scale  of  1-5,  with  1  equating  to  very  low  and  5 
equating  to  very  high  (In  terms  of  NPS  pollution  risk). 
Therefore,  substrates  within  the  watershed  corresponding  to 
H  (highly  erodlble)  are  assumed  to  be  equivalent  to  a  value 
of  5  In  terms  of  risk.   Subsequently,  substrates  that  are 
moderately  erodlble  Is  assigned  a  risk  value  of  3,  and  any 
geology  within  the  watershed  which  Is  determined  low,  with 
respect  to  erodlbillty.  Is  assigned  a  risk  value  of  1 .   In 
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other  words,  the  percentage  of  the  watershed  which  Is 
determined  as  H,  M,  or  L  Is  multiplied  by  a  risk  factor  of 
5,  3,  and  1  respectively.  These  values  are  summed  and  the 
total  for  the  watershed,  based  on  geology.  Is  determined. 
Howard  Creek  yielded  a  total  risk  value  of  1.5  and  Jones 
Meadow  yielded  a  total  value  of  2.0.  This  process  Is 
presented  In  the  formulas  found  In  Table  18  and  the  results 
are  depicted  In  Table  11. 

The  third  step  requires  risk  calculations  based  on  past 
and  future  land-use  activities.  While  this  process  may 
appear  somewhat  more  complex.  It  Is  outlined  In  the 
fol lowing  8  steps: 


Step  1:  Determine  the  year/type/acreage  of  past  and 
planned  land-use  activities. 

Step  2:  Determine  the  erosion  classification  of  the 
terrain  upon  which  the  activity  Is  to  take 
place.  This  Is  generated  by  the  GIS. 

Step  3:  Determine  slope  category  upon  which  activity 
will  or  has  taken  place  using  GIS. 

Step  4:   Obtain  risk  value  from  the  Montana 

Eros  Ion- Impact  Matrix  based  upon  slope  and 
erosion  classifications. 

Step  5:  Determine  the  percentage  of  the  watershed  upon 
which  the  activity  took  place. 

Step  6:   Obtain  a  Recovery  Coefficient  (RC)  for  the 
age/type  of  activity  using  Table  12. 

Step  7:  Multiply  the  risk  value  times  the  percentage 
of  the  watershed  times  the  RC  to  obtain  the 
total  impact  in  terms  of  risk  for  that 
activity.  (Formula  In  Table  18) 
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Step  8:  Sum  all  past  and  planned  land-use  totals. 

Once  the  total  Impact  of  land-use  activities  Is 
determined,  summation  with  risk  based  upon  slope/geology 
may  be  done.  This  Is  the  final  step  of  the  calculation 
process.  The  total  land-use  risic  values  (Table  13-B:  Howard 
Creeic  and  Table  17:  Jones  Meadow)  are  added  to  the  values 
obtained  for  slope/geology. 

The  results  for  total  risk  values.  In  terms  of 
potential  nonpoint  source  pollution,  based  on  geology/slope 
and  past  land-use  history  are  displayed  in  Table  11(A). 
Howard  Creek  Is  currently  at  a  risk  value  of  1.73  and  Jones 
Meadow,  due  mostly  to  geology.  Is  at  2.1  Therefore,  one 
might  Infer  that  Jones  Meadow  Is  at  a  higher  Inherent  risk 
than  Howard  Creek  and  would  be  more  likely  to  cause  problems 
in  terms  of  NPS  pollution  In  the  advent  of  future  land-use. 
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Land-Use  By  Year:     (Howard  Creek) 


Erosion 
Class 


Year  Of  Activity 


80       81        82       83       84        85       86       87       88 


89 


Table  4 :  Tractor  Partial  Cut  ( Hectares) 


Land-Use  By  Year:    (Howard  Creek) 


Erosion 
Class 

Year  Of  Activity 

80       81        82       83       84        85       86       87       88        89 

H 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

M 

M 

* 

* 

* 

* 

* 

* 

* 

« 

* 

L 

* 

* 

* 

« 

* 

* 

* 

* 

* 

* 

Table  5 :  Tractor  Qear  Cut  (Hectares) 


Land-Use  By  Year:    (Howard  Creek) 


Erosion 
Class 


Year  Of  Activity 


80       81        82       83       84        85       86       87       88        89 


H 


M 


^:  ■■••'' 
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Table  6 :  Cable  Partial  Cut  ( Hectares) 


Land-Use  By  Year:     (Howard  Creek) 


Erosion 
Class 

Year  Of  Activity 

80       81        82       83       84        85       86       87       88        89 

H 

* 

« 

* 

* 

* 

* 

« 

* 

* 

« 

M 

* 

« 

* 

* 

* 

« 

* 

« 
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« 

L 

* 

* 

* 

* 

* 

* 

* 

« 

* 

Table  7 :  Rubber  (Tire)  Partial  Cut  (Hectares) 


Land-Use  By  Year:     (Howard  Creek) 


EIrosion 
Class 


Year  Of  Activity 


Table  8:  Site  Prep  Dozer  ( Hectares) 


Land-Use  By  Year:     (Howard  Creek) 


80   81   82   83   84   85   86   87   88   89 


Erosion 
Class 


H 


Year  Of  Activity 


80       81        82       83       84        85       86       87       88        89 


M 


Table  9:  Broadcast  Bum  (Hectares) 


Land-Use  By  Year:     (Howard  Creek) 


Eirosion 
Class 


Year  Of  Activity 


80       81        82       83       84        85       86       87       88 


89 


H 


M 


t 


^"       'Ki 
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Table  10:  Site  Prep  LOP  ( Hectares) 


RECOVERY   COEFFICIENTS 

FOR  LAND-USE  ACIIVITIES 


LAND-USE 

YEAR  OF  ACTIVITY          | 

ACTIVITY 

89     88    87    86  85    84    83    82 

81    8ol 

TRACKED 
EOUIPMENT 

.4  5 

.4 

.35 

.3 

.25 

.2 

.15 

.1 

.05 

0 

RUBBER 
SKIDDER 

.45 

.4 

.35 

.3 

.25 

.2 

.15 

.1 

.05 

0 

PARTIAL 
SUSPENSION 

.35 

.30 

.25 

.2 

.15 

.1 

.05 

.05 

0 

0 

FULL 
SUSPENSION 

.2 

.15 

.1 

.05 

.05 

.05 

.05 

.05 

0 

0 

i 

TRACKED 
EQUIPMENT 

.35 

.30 

.25 

.2 

.15 

.1 

.05 

0 

0 

0 

RUBBER 
TIKE 

.35 

.30 

.25 

.2 

.15 

.1 

.05 

0 

0 

0 

PARTIAL 
SUSPENSION 

.2 

.15 

.1 

.05 

.05 

.05 

0 

0 

0 

0 

FUT.I. 
SUSPENSION 

.1 

.05 

.05 

.05 

.05 

.05 

0 

0 

0 

0 

1 

?9 

MACHINE 

Piip;srAR 

.8 

.8 

.7 

.6 

.5 

.4 

.3 

.2 

.1 

.05 

BROADCAST 
BURN 

.6 

.6 

.55 

.5 

.4 

.3 

.2 

.15 

.1 

.05 

OTHER 

.1 

.05 

.05 

.05 

.05 

.05 

0 

0 

0 

0 

1 
1 

PERMANENT 

1.0 

1.0 

1.0 

LO 

1.0 

1.0 

1.0 

LO 

1.0 

LO 

lEMPORARY 

.9 

.9 

.9 

.9 

.9 

.9 

.9 

.9 

.9 

.9 

Table  12 :  Recovery  Coefficients  For  Various  Land-Use  Activities 

(Klock,  1985) 


"RISK"  CALCULATIONS 
FOR  LAND -USE  HISTORY 

HOWARD  CREEK ,  LOLO  N.E 


< 

u 
>i 

x 
n 

CO 

o 

H 
EH 

< 

D 
U 

u 
\^ 

CO 
H 

^ 
g 

i 

tn 

z 

CO  CO 

o 

CO  ^ 

AREAL% 
OF  BASIN 

si 

So 

CUMULATIVE 
"RISK" 
VALUES 

1981 

SPD 

1  4 

L 

5-20 

2 

.002 

.  1 

.0004 

1983 

IPC 

14 

L 

5-20 

1 

.003 

.05 

.000  1 

5  2 

20-40 

2 

.  0  1 

.00  1 

3 

5-20 

.0006 

* 

33 

M 

5-20 

3 

.006  5 

.0009 

7  2 

20-40 

3 

.014 

.0021 

SPD 

9 

L 

0-5 

I 

.002 

.3 

.000  5 

22 

5-20 

2 

.00  4 

.00  3 

109 

M 

5-20 

3 

.02  1 

.02  5 

3  3 

20-40 

4 

.006 

.00  5 

BB 

4 

20-40 

3 

.0008 

.  2 

.0005 

1984 

TPC 

2  4 

M 

5-20 

3 

.00  4 

.  1 

.00  1 

2  3 

L 

0-5 

1 

.004 

.000  4 

180 

5-20 

2 

.03  5 

.007 

16  6 

20-40 

.032 

.006 

CPC 

33 

L 

20-40 

.0  06 

.05 

.0006 

SPD 

2  4 

M 

5-20 

3 

.00  4 

.  4 

.00  5 

2  3 

L 

0-5 

1 

.00  4 

.00  5 

17  2 

5-20 

2 

.034 

.027 

13  8 

20-40 

3 

.02  6 

.03  2 

LOP 

3  3 

20-40 

1 

.00  06 

.05 

.000  3 

Table  13  -  A :  Summary  of  Land-Use  "RISK"  Calculations 


(1980-1984) 


LEGEND 


SPD-  SITE  PREP.  DOZER 
TPC-  TRAC.  PART  CUT 
CPC-  CABLE  PART  CUT 
LOP-  SITE  PREP  LOP 


"RISK"  CALCULATIONS 
FOR  LAND -USE  HISTORY 

HOWARD  CREEK ,  LOLO  N.E 


DC 

n 
o 

H 
< 

D 

CO 
H 
0^ 

§ 
^ 

1 

(0 

o  < 

o 

(0 

AREAL  % 
OF  BASIN 

St 

§8 

CUMULATIVE 
"RISK- 
VALUES 

1985 

TPC 

3 

M 

0-5 

2 

.  0006 

.15 

.0002 

19 

5-20 

3 

.  003 

.0017 

50 

20-40 

.  0098 

.0147 

30 

L 

5-20 

2 

.  006 

.0018 

7  3 

20-40 

.014 

.0043 

SPD 

21 

L 

5-20 

2 

.004 

.  5 

.004 

32 

20-40 

3 

.006 

.009 

LOP 

31 

M 

20-40 

1 

.006 

.05 

.0003 

1986 

TPC 

5 

M 

0-5 

2 

.  0009 

.  2 

.0004 

10 

5-20 

3 

.0019 

.0011 

23 

20-40 

.  0046 

.0027 

14 

20-40 

2 

.  0027 

.001 

SPD 

5 

M 

0-5 

1 

.  009 

.  6 

.006 

10 

5-20 

3 

.  002 

.004 

24 

20-40 

4 

.  005 

.011 

10 

I. 

5-20 

2 

.  002 

.002 

4 

20-40 

3 

.  0007 

.001 

CPC 

6 

M 

20-40 

2 

.001 

.2 

0004 

ROADS 

PERMANENT 

61 

M 

20-40 

3 

.012 

1.0 

.036 

TEMPORARY 

2  1 

4 

.  004 

.012 

TOTAL 

1^ 

Table  13  -  B :  Summaiy  of  Land-Use  "RISK"  Calculations 

(1985-1989) 


Land-Use  By  Year:   (Jones  Meadow) 


Erosion 
Class 


Year  Of  Activity 


80       81        82       83       84        85       86       87       88        89 


H 


M 


Table  14 :  Rubber  Partial  Cut  ( Hectares) 


Erosion 
Class 


Year  Of  Activity 


80       81        82       83       84        85       86       87       88        89 


H 


M 


Table  15:  Tractor  Clear  Cut  (Hectares) 


Erosion 
Class 


Year  Of  Activity 


80       81        82        83       84        85       86       87       88        89 


H 


M 


Table  16 :  Site  Preparation  Dozer  (Hectares) 


"RISK"  CALCULATIONS 
FOR  LAND -USE  HISTORY 

JONES  MEADOW,  LUBRECHT  FOREST 


"RISK"    CALCULATIONS          1 
BY     YEAR                             1 

% 
^ 

1 

< 

CO 

U 

1 

Z 

o  < 

s 

CO 

Pi 

AREAL% 
OF  BASIN 

CUMULATIVE 
"RISK- 
VALUES 

1980 

SPD 

4.0 

H 

5-20 

4 

.007 

.05 

.001 

1988 

RPC 

11 

H 

5-20 

3 

.017 

.3 

.02 

SPO 

5.0 

H 

5-20 

4 

.008 

.  8 

.026 

ROADS 

PERMANENT 

4.2 

M 

5-20 

3 

.006 

1.0 

.018 

TEMPORARY 

2.1 

4 

.003 

.9 

.012 

1   TOTAL 

.08 

Table  17 :  Summary  of  Land-Use  "RISK'  Calculations 


(1980-1989) 


LEGEND 


SPD-  SITE  PREP. 

DOZER 
RPC-  RUBBER  TIRE 
PARTIAL  CUT 


"RISK"  Formulas 

H 

2 

n=L 

%WS    *   "RISK"  Value  (Geology) 

(a) 

'RISK"  Formula  Based  on  Geology 

n 

2 

i=  1 

(b)  "] 

%WS  *  "RISK"  Value  (Activity)  *  RC 

RISK"  Formula  For  Activity  Based 
on  Slope/Geology 

Table  18 :  Formulas  For  "RISK'  Calculations 
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CONCLUSIONS 

A  risk  analysis  methodology  was  developed  and  applied 
to  two  small,  third-order  watersheds  In  western  Montana. 
All  spatial  information  was  digitized  using  Geographic 
Information  Systems.  Spatial  analysis  was  conducted  using 
GIS  to  obtain  values  required  for  risk  calculation. 

Risk  values  for  Howard  Creek  and  Jones  Meadow  were 
determined  to  be  1.73  and  2.1  respectively.  These  numbers 
can  be  used  as  a  comparative  tool,  on  a  relative  basis,  to 
evaluate  watershed  sensitivity  with  respect  to  potential 
nonpoint  source  pollution  problems.  Specifically,  risk 
values  based  on  soil  erodibi I ity/slope  may  be  combined  with 
those  resulting  from  past  land-use  activity.  This  yields  a 
risk  assessment  value  for  the  entire  watershed  In  It's 
current  condition. 

After  an  initial  evaluation  Is  made  to  obtain  the 
existing  risk  value,  the  model  can  be  used  as  a  planning 
tool.  The  risk  value  may  be  Increased  by  calculating 
additional  land-use  risk  values  for  planned  activities. 

The  objective  of  the  risk  calculation  Is  to  minimize 
potential  Impacts  that  might  incur  from  planned  land 
management  activity.   Ideally,  the  overall  risk  value  for  a 
given  watershed  be, as  low  as  possible. 
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Furthermore,  areas  which  are  Innately  more  erodlble  should 
be  avoided.  The  GIS  can  produce  maps  of  sensitive  areas  and 
If  possible,  they  should  be  recognized  and  avoided  during 
planning. 

However,  It  is  left  to  the  user  of  this  methodology  to 
determine  whether  or  not  a  given  risk  value  Is  acceptable 
for  planning  purposes.  This  thesis  merely  suggests  a  means 
of  signaling  caution  for  areas  that  are  sensitive  by  nature 
or  have  been  subjected  to  extensive  land-use  In  the  past. 

The  results  of  the  initial  risk  assessment  application 
for  the  pilot  study  areas  Is  encouraging.   However, 
additional  research  is  required  to  refine  the  mechanics  of 
the  model.  Other  physical  phenomena,  such  as 
precipitation  and  delivery  distance  between  activity  and 
stream  channel,  warrant  consideration  for  possible  Inclusion 
Into  the  risk  matrix.   Furthermore,  the  development  of 
management  guidelines  concerning  resulting  risk  values  is 
necessary. 

Runoff  processes  and  cumulative  effects  should  be 
estimated  as  functions  of  soli,  vegetation,  and  land-use 
activity  and  simulated  as  past  physically  based  elements  of 
the  natural  system  on  a  temporal  and  spatial  basis  (Walsh, 
1985).  Since  watershed  monitoring  can  be  extremely 
expensive  and  time  consuming,  the  use  of  a  risk  assessment 
methodology  is  suggested. 
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Organizing  data  In  a  geographic  Information  system  and 
Integrating  risk  analysis,  permits  the  appraisal  of  nonpoint 
pollution  over  extensive  areas. 

The  capability  of  advanced  GIS  technology  allows 
editing,  data  manipulation,  and  spatial  analysis  to  be  done 
In  a  timely  and  cost  efficient  manner.  Land  management 
agencies.  Including  the  U.  S.  Forest  Service,  are  moving 
forward  In  the  adoption  of  GIS  as  a  tool  to  support  natural 
resource  management  (Bailey,  1988).   The  ability  to  turn 
spatial  Information  Into  a  digital  data  and  then  to  edit, 
store  and  display  the  data  as  maps  or  color  images  In 
combination  with  risic  analysis  makes  this  procedure  an 
invaluable  tool  for  the  management  of  water  resources. 
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APPENDIX  A;  TABLE  INTERACTION 

The  following  discussion  summarizes  the  Interaction  of 
the  Information  contained  In  Tables  2-17.  Tables  2  and  3 
are  used  to  calculate  the  risk  value  of  the  watershed  based 
so  ley  on  geology.  The  results  are  depicted  In  Table  11. 
Tables  4-16  Include  a  summary  of  land-use  activity,  by  year, 
for  both  pilot  watersheds.  The  land-use  values  are 
multiplied  by  recovery  coefficients  In  Table  12,  and  the 
results  are  summed  In  Table  13/13A  and  17  for  Howard  Creek 
and  Jones  Meadow  respectively.  Total  risk  values  for 
geology  and  land-use  are  combined  and  the  resulting 
"cumulative"  risk  value  for  both  watersheds  are  shown  In 
Table  11A  (See  Figure  14). 
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A  FOREST  MANAGEMENT  NONPOINT  SOURCE  RISK  ASSESSMENT 
GEOGRAPHIC  INFORMATION  SYSTEMS  APPLICATION 
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ABSTRACT:  The  map-based  nonpoint  source  risk  assessment  methodology 
developed  by  the  State  of  Oregon  for  their  208  Assessment  has  been 
adapted  for  use  with  Geographic  Information  Systems.  Central  to  the 
CIS  Application  is  a  "Leopold-type"  risk  matrix  that  has  been  developed 
concurrently.  The  orthogonal  matrix  displays  combinations  of  forest 
management  practices  and  classes  of  slope  and  soil  erodibility.  The 
relative  nonpoint  source  risk  for  each  matrix  cell  was  assessed  by 
experts  using  the  Nominal  Group  Technique.  The  application  of  the 
procedure  for  cumulative  effects  analysis  in  Howard  Creek  on  the  Lolo 
National  Forest  is  presented. 

KEY  TERMS:  Cumulative  watershed  effects;  risk  assessment;  nonpoint 
source;  forest  management. 

INTRODUCTION 

The  assessment  of  nonpoint  source  pollution  and  cumulative  watershed 
effects  resulting  from  forest  management  is  proving  to  be  a  difficult 
task.  Application  of  Best  Management  Practices  (BMPs)  does  not  provide 
an  assessment  of  condition,  monitoring  can  be  prohibitively  expensive, 
and  our  quantitative  modeling  skills  are  in  their  infancy.  An 
alternative  approach,  watershed  risk  analysis,  may,  however,  provide 
resource  managers  with  a  powerful  tool  to  assist  in  land-use  decision 
making.  In  addition,  Geographic  Information  Systems  (GIS)  allow  the 
correlation  of  land  cover  and  topographic  information  such  as  terrain 
configuration  and  drainage  networks,  thus  making  GIS  useful  in  assessing 
the  potential  effects  of  land-use  activities  on  water  resources  (Walsh, 
1985) . 

We  have  developed  a  methodology  to  assist  in  assessing  forest 
management-related  nonpoint  source  pollution  and  cumulative  watershed 
effects  that  combines  watershed  risk  assessment  with  GIS.  Having  a  GIS 
make  a  map  was  the  easy  part  of  the  process.  The  difficult  part  was 
telling  the  GIS  which  landscape  attributes  were  to  be  on  the  map.  The 
key  to  the  application  was  the  construction  of  a  nonpoint  source 
pollution  risk  matrix.  The  results  of  the  application  on  a  heavily 
impacted,  mixed-ownership  watershed  are  presented. 
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RISK  MATRIX  CONSTRUCTION 

Leopold  and  others  (1971)  were  among  the  first  to  propose  the  use  of 
an  environmental  matrix  to  evaluate  potential  environmental  impact. 
The  matrix  provides  a  comparatively  simple  system  intended  to  be  a  guide 
in  the  impact  assessment  process.  This  approach  was  adopted  by  Rickert 
and  others  (1978)  to  evaluate  erosion  potential  for  the  Oregon  208 
nonpoint  source  assessment  project  and  by  Brown  III  and  others  (1979) 
as  a  guide  for  land-use  planning. 

The  orthogonal  "risk"  matrices  developed  in  these  studies  assessed 
the  relative  impacts  of  various  land-use  activities  on  different 
combinations  of  slope  and  substrate  erodibility.  Forest  management 
practices  were  among  the  land-use  activities  considered,  but  Brown  III 
and  others  suggested  that  matrices  should  be  developed  regionally  to 
reflect  local  conditions,  practices,  and  socio-economic  considerations. 

Therefore,  we  contacted  a  large  group  of  soil  and  water  specialists 
and  silviculturists,  representing  land  management  agencies  and  the 
forest  products  industry  in  Montana,  to  build  the  Montana  erosion- 
impact  matrix.  We  quickly  learned  that  it  was  not  going  to  be  as  easy 
as  we  thought.  There  was  considerable  disagreement  over  the  relative 
risk  values,  on  a  simple  scale  of  1  to  5,  to  be  used  in  the  body  of  the 
matrix  as  well  as  the  basic  combinations  of  slope  and  erodibility  that 
would  produce  different  potential  impacts. 

We  convened  the  group  of  experts  twice,  and  through  the  Nominal  Group 
Technique,  arrived  at  the  consensus  matrix  appearing  in  Figure  1.  The 
horizontal  axis  is  composed  of  4  principal  forest  management  practices: 
1)  clear  cutting;  2)  partial  cutting;  3)  site  preparation;  and  4)  road 
construction.  Each  principal  practice  is  further  split  into  specific 
methods,  treatments  or  design. 

The  group  decided  that  4  slope  classes  on  the  vertical  axis 
were  sufficient  ,and  distinctly  different  in  inherent  erosion  risk,  and 
that  3  basic  soil  erodibility  classes  were  also  adequate.  These  classes 
were  based  primarily  on  the  geologic  erosion  factors  described  in  the 
R1-R4  Sediment  Yield  Prediction  Procedure  (USDA  Forest  Service  1981) . 
Everyone  agreed  that  alluvial  and  granitic  soils  were  highly  erodible, 
those  formed  from  Belt  Series  metamorphics  had  low  erodibility,  and  that 
most  other  substrates  fell  into  the  moderate  erodibility  class. 

In  all  but  19  of  the  156  cells  in  the  matrix,  the  group  accepted  the 
modal  response  for  consensus.  Those  cells  marked  with  an  asterisk  still 
have  to  be  resolved,  and  this  produces  some  obvious  discrepancies  in  the 
matrix.  For  example,  partial  cutting  on  highly  erodible  0-5%  slopes  has 
more  risk  than  clear  cutting.  Regardless,  the  completed  matrix  provided 
the  various  attributes  to  be  mapped  by  the  CIS  for  this  application  - 
4  slope  classes,  geologic  substrates,  areas  clear  cut  or  partial  cut 
by  4  different  methods,  areas  receiving  site  preparation,  and  roads. 

HOWARD  CREEK 

The  Howard  Creek  watershed,  totalling  5015  hectares,  is  located 
approximately  32  kilometers  southwest  of  Missoula,  Montana.  Elevation 
in  the  watershed  ranges  from  1190  to  1770  meters.  Average  annual 
precipitation  at  mean  elevation  is  about  100  cm,  and  roughly  60%  of  that 
falls  as  snow.  The  forest  cover  is  primarily  Douglas-fir  mixed  with 
sub-alpine  fir  at  the  higher  elevations  and  Ponderosa  pine  at  the  lower 
elevations  and  south  aspects.  The  main  stem  of  Howard  Creek  is  oriented 
east-west  while  the  3  sub-basins.  Tepee  Creek,  North  Fork  of  Howard 
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Creek,  and  Krystal  Creek,  are  oriented  north-south  (See  Figure  2).  This 
map,  and  the  following  maps  were  prepared  using  PAMAP  (1989)  CIS 
software  and  information  hand-digitized  at  a  scale  of  1:24000. 

The  Lolo  National  Forest  Land  Systems  Inventory  maps  the  primary 
geologic  groups  found  in  Howard  Creek  as  metasedimentary  and 
undifferentiated  materials  (see  Figure  3)  . 

Metasedimentary  parent  materials  are  derived  from  Belt  Super  Group 
quartzite,  argillite,  and  siltite.  Rock  fragment  hardness  is  variable 
depending  upon  the  degree  of  rock  weathering.   Weathering  is  dependent 
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Figure  1.   The  Erosion-Impact  Risk  Matrix. 
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Figure  2.   The  Howard  Creek  Watershed,  Montana, 


Figure  3.   Geology  in  the  Howard  Creek  Watershed. 


26 


on  associated  faults,  preponderance  of  argillites,  and  calcium  carbonate 
content.  These  materials  were  classified  as  either  L-low  or  M-moderate 
erodibility  for  the  erosion-impact  matrix. 

Undifferentiated  geology  is  composed  of  materials  derived  from  Belt 
Super  Group  metasedimentary  rocks  or  weakly  weathered  granitic  rocks. 
Materials  include  alluvium  on  terraces  and  flood  plains;  shallow  soils 
on  flood  scoured  foot  slopes  and  stream  breaklands,  strongly  frost 
churned  broadly  convex  ridges,  and  glacial  outwash  on  plains  These  are 
classified  as  H-highly  erodible  in  the  Montana  erosion-impact  matrix. 

Like  much  of  western  Montana,  the  watershed  has  a  "checkerboard" 
pattern  of  land  ownership.  Champion  Timber lands  owns  34%  of  the 
watershed.  Plum  Creek  Timber  owns  23%  of  the  watershed,  and  the 
remaining  42%  is  managed  by  the  Lolo  National  Forest.  Ownership  was 
mapped  by  the  CIS  and  is  shown  in  Figure  4.  Timber  harvest-related 
activities  between  1981  and  1986  impacted  17%  of  the  watershed.  The 
location  of  these  activities  is  shown  in  Figure  5.  In  a  report  to  the 
Lolo  Forest  Supervisor  (Munther  and  others  1987) ,  cumulative  watershed 
effects  in  Howard  Creek  were  estimated  to  have  produced  a  sediment  load 
increase  of  50%  and  water  yield  increase  of  8%.  The  Forest  Supervisor 
immediately  imposed  a  10-year  timber  harvest  moratorium  on  Forest 
Service  lands.  Industry  has  cooperated,  and  there  has  been  little 
activity  in  the  watershed  during  the  past  3  years. 

A  CUMULATIVE  EFFECTS  RISK  INDEX  VALUE 

An  undisturbed  watershed  has  a  0  risk  index  value.  The  risk  index 
value  for  a  managed  watershed  depends  on  the  areal  extent  of  activities, 
the  age  of  the  treatments,  and  the  type  of  terrain  on  which  they  are 
located.  The  calculation  of  a  cumulative  effects  risk  index  value  is 
straightforward  and  is  outlined  in  the  following  8  steps: 

Step  1:   Determine  the  year/type/acreage  of  past  (or 

planned)  forest  management  activities.   Permanent 
road  disturbance  is  calculated  at  .87  ha/km; 
temporary  road  disturbance  is  calculated  at 
.37  ha/km.   Digitize  this  information  for  use  in 
the  CIS  if  it  has  not  already  been  done. 

Step  2:   Determine  the  erosion  classification  of  the 

terrain  on  which  the  activity  has  taken  (or  will 
take)  place.   This  is  generated  by  the  CIS. 

Step  3:   Determine  slope  category  on  which  the  activity 
has  taken  (or  will  take)  place  using  CIS. 

Step  4:   Obtain  a  risk  value  for  each  activity  from  the 
Erosion-Impact  Matrix  that  has  been  developed. 

Step  5:   Determine  the  areal  percentage  of  the  watershed 
on  which  the  activity  took  place. 

Step  6:   Use  a  recovery  coefficient  (RC  -  a  value  between  1 
and  0)  to  reduce  the  risk  associated  with  older 
activities.   In  this  application,  we  assumed  a 
linear  10-year  risk  recovery  for  all  activities 
except  roads,  which  always  maintain  a  RC  of  1. 
Thus,  current  year  activities  have  a  RC  of  1, 
4  year-old  activities  have  a  RC  of  0.6,  and  so  on. 

Step  7:   Multiply  the  risk  value  times  the  percentage 
of  the  watershed  times  the  RC  to  obtain  the 
total  risk  for  that  activity. 

Step  8:   Sum  all  past  and  planned  land-use  totals. 
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Figure  4.   Checkerboard  Ownership  in  the  Howard  Creek  Watershed 


Figure  5.   Timber  Harvest  Activity  in  the  Howard  Creek  Watershed. 
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Table  1  contains  a  summary  of  disturbance  information,  by  year,  for 
Howard  Creek.  Harvest,  site  preparation  and  road  data  were  obtained 
directly  from  the  3  land  owners  in  the  watershed.  To  avoid  double 
accounting,  only  the  last  entry  onto  a  site  was  counted.  Usually  this 
was  for  site  preparation,  which  was  almost  always  dozer  piling  of  slash. 
If  no  site  preparation  had  been  completed,  then  harvesting  was  the  only 
disturbance  counted,  and  this  was  almost  always  partial  cutting  with 
tractor  skidding. 

Table  1.   "Risk"  Calculations  for  Howard  Creek,  Lolo 

National  Forest,  Montana. 


%  of 

1990 

1986 

Area 

total 

weighted 

weighted 

Year 

Activity 

(ha) 

area 

risk 

risk 

1981 

Site  Prep. 

14 

0.27 

.0004 

.0020 

1983 

Site  Prepo 

177 

3.53 

.0342 

.0798 

1984 

Harvest  only 

36 

0.72 

.0065 

.0130 

Site  Prep. 

390 

7.78 

.  .0650 

.1300 

1985 

Harvest  only 

92 

1.83 

.0275 

.0495 

Site  Prep. 

84 

1.67 

.0160 

.0288 

1986 

Harvest  only 

6 

0.12 

.0021 

.0035 

Site  Prep. 

53 

1.06 

.0247 

.0412 

Permanent  Roads 

61, 

,4 

1.22 

.0360 

.0360 

Temporary  Roads 

21. 

.3 

0.42 

.0160 

.0160 

Unrecovered  Disturbance 

18.62% 

1990  Cumulative 

"Risk" 

.2284 

1986 

Cumulat: 

Lve  "Risk" 

.3998 

Note  that  in  1990,  the  cumulative  risk  index  value  is  0.2284,  but  if 
the  risk  index  value  had  been  calculated  in  1986,  the  year  the  timber 
harvest  moratorium  was  called,  it  would  have  been  0.3998.  The 
differences  between  the  1990  and  1986  weighted  risk  columns  in  Table  1 
are  due  only  to  the  effect  of  the  recovery  coefficients. 

THRESHOLDS? 

How  much  disturbance  a  watershed  can  tolerate  before  reaching  a  level 
that  produces  significant  environmental  damage  is  both  technically  and 
politically  difficult  to  answer.  In  the  Forest  Service's  Region  5,  a 
watershed's  "Threshold  for  Concern"  is  typically  an  areal  disturbance 
(actually,  equivalent  road  acres)  of  between  10  and  20%  depending  on  the 
watershed's  sensitivity  (Cobourn  1989).  These  values  are  based  on 
research  conducted  by  Harr  and  others  (1975)  . 

If  Howard  Creek  were  in  Region  5,  its  combination  of  slopes  and 
geology  would  probably  classify  it  as  moderately  sensitive,  and 
therefore  would  have  a  "Threshold  for  Concern"  of  15%.  Note  in  Table 
1  that  6  years  of  management  produced  an  unrecovered  disturbance  of 
nearly  19%.  The  "yellow  flag  of  caution"  would  not  yet  have  been  raised 
because  of  recovery  and  only  some  percentage  of  the  disturbance  was 
equivalent  road  acres.  The  specialists  on  the  Lolo  National  Forest 
raised  their  own  "yellow  flag  of  caution"  based  on  monitoring 
information.  Nevertheless,  there  were  some  indications  of  changes  in 
the  system,  so  perhaps  a  10  to  15%  disturbance  threshold  is  legitimate. 
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'  Similar  rational  for  assigning  a  risk  index  value  threshold  may  be 
used.  A  moderate  risk  for  generating  cumulative  effects  has  a  value  of 
3  in  the  Erosion-Impact  matrix.  If  we  were  to  chose  a  15%  areal 
disturbance  limit,  and  be  willing  to  accept  moderate  risk,  then  the  risk 
index  target  would  be  0.45  (0.15  X  3) .  A  10%  disturbance  limit  would 
have  a  risk  index  of  0.3.  Similarly,  acceptance  of  higher  or  lower  risk 
could  also  increase  or  decrease  the  target.  In  either  event,  the  Howard 
Creek  watershed  was  being  stressed,  or  nearly  so,  in  1986.  Recovery  is 
rapid,  however,  so  deferring  activity  or  selecting  low  risk  options  (a 
combination  of  sites  and  methods)  would  allow  the  cumulative  risk  value 
to  fall.  Note  that  by  1990,  the  cumulative  risk  index  value  fell  to 
about  0.23. 

FUTURE  EFFORTS 

We  are  about  to  begin  testing  this  risk  assessment  methodology,  using 
the  GIS  in  operation  on  the  Flathead  National  Forest,  as  part  of  a 
comprehensive  study  to  assess  cumulative  watershed  effects  in  the 
Flathead  Basin  funded  by  the  Flathead  Basin  Fisheries  and  Water  Quality 
Cooperative.  Additional  projects  are  already  underway  to  use  the 
nominal  group  technique  to  try  to  obtain  expert  consensus  on  both 
recovery  rates  and  disturbance  thresholds. 
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MONTANA  "RISK"  MATRIX  (3/9i) 


TYPE 


CLEARCUT 


PARTIAL  CUT 


SITE  PREP 


ROADS 


SLOPE 


SOIL 
EROD. 


^ 


(X4 


"^5 


pCQ 

0413 

■<  m 


E- 

OS  pa 
n 


K 

Eh 

o 


a, 


H 


0-5% 


M 


>** 


)** 


)** 


I** 


2** 


>** 


1     1 


H 


5-20% 


M 


4* 


2* 


H 


20-40% 


M 


5* 


I** 


1*4 


4* 


3* 


H 


>  40% 


M 


(  *  Indicates  Mean  Value  Used,  **  Final  Group  Consensus  Adyustment  for  Internal  Ck>n8istency) 
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cost  of  $7.00  per  copy,  for  a  total  cost  of  $350.00,  which  includes 
$350.00  for  printing  and  $0.00  for  distribution. 
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